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TECHNICAL NOTE 3830 

GROWTH OF DISTURBANCES IN A FLAME -GENERATED SHEAR REGIONI 
By Perry L. Blackshear^ Jr. 


SUMMARY 

The growth of transverse velocity disturbances in the shear region 
caused by a flame in a duct is examined theoretically and experimentally. 

The theoretical analysis consists of a stability analysis of a flow 
field arising from a flame in a duct. The flow is found to be neutrally 
stable to symmetric disturbances and unstable to antisymmetric ones. 

In following the history of a disturbance of given frequency from 
the flameholder downstream^ it is found that a most -amplified frequency 
exists, which is equal to flow velocity divided by duct width. Disturbances 
at frequencies greater or less than this are amplified less than this one. 

The amount that an incoming disturbance is amplified depends on the 
flame slope. For the most -amplified frequency, the ratio of terminal to 

initial transverse disturbance velocity is equal to where 

the slope is roughly turbulent flame speed divided by flow velocity. 

This ratio of terminal to initial transverse disturbance velocity is de- 
pendent on inlet velocity and flame speed only in the stated manner and 
is essentially independent of density ratio. 

The critical wavelength separating disturbances that grow from those 
that do not grow varies as a constant times flame width, where the con- 
stant is 2“ for flames that are narrow compared with duct width and be- 
comes infinite as the flame reaches the wall. The flow is unstable to 
disturbances having wavelengths greater than critical. Thus, short- 
wavelength (i.e., high-frequency) disturbances achieve their terminal 
amplifications near the flameholder, while long -wavelength disturbances 
achieve theirs farther downstream. The most -amplified disturbance 
achieves its terminal amplification when the flame fills 0-60 of the duct. 

^Originally submitted as a doctoral thesis at Case Institute of 
Technology under professor Harold G. Elrod. 
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The results are related to measurable quantities in the form of 

(1) Flame-front displacement 

(2) Phase velocity 

(3) To a first approximation, the effect of disturbance frequency 
on flame propagation velocity 

The experimental part of the program consisted of imposing disturb- 
ances of various frequencies upon a flame stabilized in a duct and measur- 
ing the effects through shadow photography and photomultiplier probe sur- 
veys. The efforts to verify elements of the theory are encouraging but 
not conclusive. Reasonable agreement is achieved in the comparison of 
flame-front displacement and phase velocity. In the measurements of the 
effect of frequency of disturbance on flame speed, there appeared to be 
a most -effective frequency, just twice the most -amplified frequency given 
by theory. 

There is no direct evidence of transition due to those disturbances. 
The flow was turbulent in all cases owing to transition in the flameholder 
boundary layer, either on the flameholder or near the point of separation. 
There was evidence, however, that the base turbulent flame, wrinkled by 
the long-wavelength disturbances, was altered (made more turbulent) when 
the disturbances were present. In other words, the large-amplitude, 
long-wavelength disturbances cause a sort of transition in the already 
turbulent flow, amplifying the smaller scale disturbances present. 


INTRODUCTION 

This work investigates the stability of the flow field arising from 
an anchored flame in a duct, a subject which is of interest for two 
reasons: 

(1) Experimental investigations have shown that the turbulent flame 
speeds attained in a flame anchored in a duct depend more strongly on the 
turbulence originating in the shear regions generated by the flameholder 
and flame than on the turbulence in the incident gas stream. 

(2) Some transverse modes of burner resonance encountered have been 
accompanied by an increase in the rate of flame propagation. 


DEPENDENCE OF TURBULENT FLAME SPEED ON DISTURBANCES 
ORIGINATING IN SHEAR REGIONS 

The problem of describing the flow field that arises when a flame is 
anchored in a duct was first posed and solved by Scurlock (ref. l). One 
result of this work was the relation between the fraction of mixture 
burned and the local flame width. Scurlock used this result to measure 
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turbulent flame speeds in an experiment originally designed to give the 
effect of inlet turbulence on flame speeds. He found that turbulence 
associated with two regions of shear flow masked, to a great extent, the 
effect' of imposed stream turbulence. These two regions consist of the 
near wake of the flameholder and the shear region associated with the 
fully developed flame in a duct. The latter shear region is caused by 
the pressure gradient acting on gases of different densities. Scurlock 
calculated the velocity profile in this region. 

In flow past cylinders without flames, eddies are encountered. 
Scurlock found that the flame seated on a cylindrical flameholder pre- 
vented eddies from shedding from the flameholder over the range of 
Reynolds numbers investigated (up to approximately 15,000). Haddock 
(ref. 2) and Zukowski (ref. 3) in investigating the mechanism of flame 
seating on cylindrical rods found that the boundary layer separating 
from either side of the rod underwent transition from laminar to turbu- 
lent flow. The transition point was found to approach the point of sep- 
aration at a Reynolds number of about 10^. The work of references 2 and 
3 shows that transition in the flameholder boundary layer is one mechanism 
for producing turbulence in the shear region near the flameholder. 

The role played by the second shear region cited by Scurlock as a 
source of turbulence has not been directly investigated. After obtain- 
ing experimental results similar to those of Scurlock, several investi- 
gators (e.g., refs. 4 and 5) have arrived at the conclusion that turbu- 
lence is generated in this region. The types of argument leading to 
these conclusions are briefly indicated in the following discussion. 

Theories of turbulent flame propagation have been postulated for . 
two types of turbulence: 

(1) Turbulence with a scale small compared to a laminar flame-front 
thickness 

(2) Turbulence -with a scale large compared with laminar flame-front 
thickness (refs. 6 to 8) 

Experiments conducted in the absence of shear regions gave good agreement 
between values predicted by the large-scale turbulence theory and experi- 
mental values of flame speed (refs. 9 to 12). Experiments conducted in 
the presence of shear regions (refs. 5, 8, and 13) gave flame speeds 
higher than those predicted by theory. 

There is a current controversy concerning the interpretation of 
this difference (refs. 5, 8, and 13 to 16). There is general agreement 
that, when a shear region exists and the shear region in question is of 
the type associated with an anchored flame in a duct, the higher flame 
speed is due in part to turbulence amplified by or originating in this 
shear region. 
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Rayleigh (ref. 17) observed that some jets and diffusion flames are 
sensitive to certain imposed transverse disturbances. He performed an 
inviscid stability analysis on a -_y\_ velocity profile and found the 
profile unstable to antisymmetric disturbances, but neutrally stable to 
symmetric wavy disturbances. This result agreed with his observations 
on sensitive flames and jets. Some of the limitations of Rayleigh's • 
work are considered in the section THEORY; STABILITY OF FLOW FOR A FLAME 
IN A DUCT . The analytical method used by Rayleigh is applied in the 
present examination of the instability of flame-zone flow to imposed 
disturbances . 


EFFECT OF COMBUSTOR OSCILLATIONS IN INCREASING 
RATE OF FLAME PROPAGATION 

A great variety of combustion-chamber resonances have been observed 
in jet-engine combustors (refs. 18 to 21). These resonances most often 
involve some mode of acoustic oscillation in the combustion chamber. 
Reference 18 discusses modes that involve longitudinal gas and flame 
motion and reports a falloff in combustor performance when these oscilla- 
tions occur. 

Reference 19 discusses a type of resonance identified as a transverse 
mode of oscillation in the combustion chamber. Accompanying this type of 
oscillation was an apparent increase in combustion efficiency (ref. 19) 
and in the rate of flame spreading (ref. 20). There is a possibility 
that this type of oscillation with an acoustic wavelength controlled by 
duct dimensions could excite a wavy disturbance that would subsequently 
be amplified by the flame -generated shear region. 

Unfortunately, these transverse modes (called screech) also cause 
the destruction of the burner walls when the screech becomes severe 
enough (ref. 20). In reference 21, pictures are shown of the flame near 
the flameholder when a combustor was resonating in a' transverse mode. A 
strong interaction was apparent in the form of flame wrinkles forming at 
the flameholder and growing in amplitude as they moved downstream. It 
seems apparent that at least one of the ways the transverse oscillation 
increases the rate of flame spreading is by wrinkling, hence extending, 
the interface between the burned and unburned gases. 

If the oscillation that disturbs the flame at the flameholder is 
supplied locally (by some suitable exciting mechanism) at the flameholder 
without relying on the sometimes destructive combustor resonance to do 
the exciting, some of the benefits of screech can be realized without its 
penalties. The gains in flame speed should be greatest for those fre- 
quencies of excitation that are amplified by the flame. It is hoped that 
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by examining the sensitivity of a flame in a duct to imposed disturbances 
an excitation frequency may be found that gives the optimum increase in 
flame speed. 

This report is divided into three parts, theory, experiment, and 
discussion. A separate introduction is given for the first two parts 
describing the intent. The discussion weighs the results of the other 
parts and contains some extrapolations. 


SYMBOLS 


The principal symbols are defined in the following list; special 
symbols are defined where used. 


(AdAu)a 


ratio of turbulent undisturbed flame siLrface area to disturbed 
flame surface area measured from flame photographs 


(ApAu)? 


ratio of turbulent undisturbed flame surface area to disturbed 
flame surface area measured by photomultiplier 


a. A, B, 
C,D 

b 

b' 

c 

‘^i 

Cr 

F 

f 

g 

g’ 

h 

^1 


constants of integration 

width of sheax region in Rayleigh's problem 
half -width of plateau in velocity profile 
complex phase velocity, ft/ sec • 
imaginary part of complex phase velocity 

real part of complex phase velocity 
fraction burned 
frequency, cps 

1 

gain; v» = v^ g 
gain; v' = v^ g' 

flame-front displacement from mean, ft 
duct half -width, ft 



Ki 


Kc 


Lx 

If 

n 

m 

P 

P 

P* 

Re 

Sl 

t 

u 

Uq 

Ul 

u' 


MCA TN 3830 


flame half -width, ft 
imaginary number, 
dli2 

constant, for an assumed flame-spreading mechanism 

constant, 

' V 


Eulerian scale of turbulence, ft 
length of flame measured from wake, ft 


len^h of flameholder wake, ft 

rj + rl 

- t2 

mean pressure, Ib/sq ft 
pressure, Ib/sq ft 
perturbation pressure, Ib/sq ft 
Reynolds number 
laminar flame speed, ft/sec 

turbulent flame speed, ft/sec 
time, sec 

flow velocity in x-direction, ft/sec 

flow velocity in x-direction at plane of flameholder, ft/sec 

flow velocity in x-direction in cold gas, ft/sec 

flow velocity in x-direction at apex of hot-gas velocity profile, 
ft/sec 

disturbance velocity in x-direction, ft/sec 
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v’ disturbance velocity in y-direction, ft/sec 

X coordinate in direction of duct axis (direction of mean flow), 

y coordinate at right angle to direction of mean flow, ft 

a wave number of disturbance, ft"^ 

Oq wave number of disturbance based on initial velocity and 

frequency, f2rt/UQ 


ag wave number for neutral stability based on inviscid analysis, 
p hot-gas velocity gradient, (U* - U]^)/h 2 , sec"^ 


-oh-| 

e 


-oho 

e 




1 + ri 



X 


P 

P* 

Pi 


<p(y) 


wavelength^ ft 
— , ft 

density, lb sec^/ft^ 
perturbation density, lb sec^/ft^ 
density of cold gas, lb sec^/ft^ 

density of hot gas, lb sec^/ft^ 

y-distribution of disturbance amplitude, ft^/sec 
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THEORY: STABILITY OF FLOW FOR A FLAME IN A DUCT 

INTRODUCTION 

The Steady-State Flow Field of a Flame Anchored in a Duct 

The problem of describing the flow field that exists when a flame 
is anchored in a duct was first posed and solved by Scurlock in 1948 
(ref. l). Subsequently, Tsien (ref. 22) and others (ref. 23) have elab- 
orated on the problem. The following table summarizes Scurlock' s and 
Tsien* s assumptions and results as they pertain to the present work. 


Assiimptions 

Scurlock 

Common 

Tsien 

1. Fluid is incompres- 
sible. 

1. Fluid is steady- 
state, two-dimensional, 
inviscid, non-heat- 
conducting, non- 
turbulent . 

1. Two cases are treat- 
ed, the compressible 
and incompressible. In 
the incompressible case 
Scurlock* s assumption 
(2) holds; in the com- 
pressible case, density 
is allowed to vary in 
the direction of flow. 

2. Density jumps from 
Pi to p 2 across the 
flame front, is con- 
stant in each region. 

2. Velocity, density, 
and pressure are uni- 
form at combustion 
inlet . 

2. The density ratio is 
a constant across the 
flame front at any 
station. 


3, Static pressure is 
uniform at any cross 
section. 

3. The velocity profile 
is composed of four 
straight-line segments. 


4. The flameholder is 
of negligible size, 
spans the combustor. 



5. The flame front is 
of negligible thickness. 



6. The chamber is of 
constant width. 


Results 

1. Velocity profile as 
a function of flame 
width: (slightly 

concave to incoming 
gas. ) 

1. Fraction burned as 
a function of flame 
width 

1. Peak hot-gas veloc- 
ity as a function of 
flame width 

2. Momentum pressure 
drop as a function of 
flame width (i.e., of 
fraction burned) 

2. Cold-gas velocity 
as a function of flame 
width 

2. Limiting inlet ve- 
locity beyond which 
complete combustion 
cannot occur as a func- 
tion of density ratio 
(for compressible case) 
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A stability analysis is made of a parallel flow having a velocity 
profile approximating the one found by Scurlock by means of straight- 
line segments and a density distribution as taken by Scurlock, that is, 
constant 9 \/^ 2 * 


The Stability Problem 

The study of the origins of turbulent motion has received consider- 
able attention in the past, particularly with respect to turbulence oc- 
curring in boundary layers. One of the ways the problem has been treated 
is to assume the existence of infinitesimal wavy disturbances in the 
flowing stream and seek the conditions under which these disturbances are 
amplified (refs. 24 and 25). In the application of this method to the 
study of the stability of the laminar boundary layer, a major objective 
has been to determine a minimum Reynolds number below which all disturb- 
ances decay and the effect of velocity profile on this minimum Reynolds 
number. 

The relation between the amplification of wavy disturbances and the 
transition to turbulence is supposedly that one leads to the other under 
certain circumstances. The desire to maintain laminar flow on airfoils 
makes an understanding of transition an important practical problem. In 
1941, Schubauer and Skramsted (ref. 26) conducted an experiment in which 
the disturbances introduced into a laminar boundary layer were carefully 
controlled. They found a relation between the wave numbers of amplified, 
neutral, and damped distiirbances and the Reynolds number that confirmed 
the predictions of the infinitesimal wavy disturbance theory in almost 
every important respect. 

The practical value of learning about the stability of a flow field 
created by a flame in a duct would be in fiirnishing information that 
would enable combustor designers to encourage the growth of disturbances 
in this region, because turbulence generally aids flame spreading. The 
type of information needed is (l) a range of Reynolds numbers over which 
the flow is unstable and (2) the type of disturbance that possesses the 
greatest rate of amplification within this range. In the following 
qualitative discussion it appears that most flames of practical interest 
are in a Reynolds nmber range (l) where amplified disturbances always 
exist and (2) where wavelengths of disturbances that will be amplified 
may be obtained, neglecting viscosity. 

If the disturbance varies as 

v‘ = (p(y)ei°'(^-‘=^) 

where the wave number a is real^ and the phase velocity c is gsner- 
ally complex, then, as shown by Lin (ref. 24), the stability of flows 
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of the boundary-layer type (or of flows with symmetric profiles in a two- 
dimensional channel) can be illustrated by a graph of values of a for 
neutral stability (i.e.^ where the imaginary part of c vanishes) against 
Reynolds number Re. Two types of profile are discussed. 

A profile without a point of inflection would possess the following 
plot of a against Re for Cj^ = 0: 



This curve is called the neutral envelope for the profile in question. 
Inside the envelope^ disturbances are amplified (cj^ > O) ; outside^ dis- 
turbances are damped (c^^ < O). Proceeding to the right from the minimum 
Reynolds number attained by the Cj^ = 0 curve are two branches of the 
curve Cj^ = 0. As Re becomes indefinitely great, both asymptotically 
approach a = 0. 

For the profile with a point of inflection, a similar' relation 
a(c,Re) is obtained, as shown in the following graph 
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The significant difference hetween this plot for a profile with an 
inflection point and the preceding plot for a profile without an inflec- 
tion point is that the asymptotic values of the two branches approach 
different limits for large values of Re. The upper branch approaches 
a = the wave number for neutral stability obtained by neglecting 

viscosity; the lower branch approaches a = 0 as it does for a profile 
without an inflection point. The Reynolds number range over which the 
flow with an inflection point is unstable extends to indefinitely large 
Reynolds numbers. The flow generated by a flame in a duct has, in effect, 
an inflection point. 

It is desirable either by example or direct calculation to determine 
the order of magnitude of the Reynolds number range over which the neutral 
a departs from ag, that is, where viscosity has a dominating influence. 

In two examples in the literature, Pretsch (ref. 27) and Lessen (ref. 28) 
examined a flow with a flex for stability. In both instances the minimum 
Reynolds number for Cj_ = 0 was of the order of 10; the maximum Re, 

where a for Cj_ = 0 differed from ag calculated from inviscid theory 
by more than 2 percent, was of the order of 10^. The profile of refer- 
ence 27 is for a boundary layer with. pressure rise; that of reference 28 
for the free boundary layer between streams of differing velocities. 

If these values are taken to suggest the order of Reynolds number 
where the stability of a profile with a flex is sensitive to the effects 
of viscosity, some idea of flame-zone sizes can be estimated for the 
range of velocities encountered in jet engines, that is, 50 to 500 feet 
per second. For example, at 50 feet per second, the flame width where 
viscous forces would be important would range from 1/3000 to l/30 inch. 

At 500 feet per second, the range is l/30,000 to l/300 inch. Most of the 
flame zones of practical interest are much wider than these dimensions 
and may be considered well within the range of Reynolds numbers where 
disturbances are amplified and where viscosity effects are small. There- 
fore, there is reason to expect the following inviscid stability analysis 
to give meaningful information on the value of the wave numbers a where 
amplification can occur. 

Lin (ref. 24) and Heisenberg (ref. 29) have criticized the use of 
broken-line profiles because there is a separate solution for c for 
each corner of the profile. These objections might be met by the follow- 
ing observations: 

(l) The steady-state flow field solved by Scurlock possesses corners 
(places where the velocity gradient jumps) and can be very closely approx- 
imated by straight-line segments without introducing new corners into the 
profile. The jumps in velocity gradient that occur at the corners result 
from the assumption that the flame front is of negligible thickness com- 
pared to the width of the flame zone. 
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(2) In the results that follow it is found that the flame zone is 
unstable to disturbances having wavelengths 2.5 times the flame-zone 
width and greater. Thus, any rounding of the corners over a distance 
that is small compared to the flame -zone width should not affect the re- 
sults; to the long wavelength disturbances that are amplified, the cor- 
ners will still appear sharp. 

To examine by example the effect of corners in a velocity profile, 
three profiles of a similar shape are shown in the following sketch. 

Each of these profiles has been examined for instability and a wavelength 
for neutral stability has been computed. The scale of the abscissa is 
chosen so that the value of the neutral wavelength is the same for the 
three profiles. 



The two broken-line profiles are from Rayleigh (ref. 17); the continuous 
profile from Savic (ref. 30) . These solutions correspond to an anti- 
symmetric disturbance. The flow in each case is in viscid. 


A similar comparison can be made between the velocity profile 


U = 1 + sin 


- 1 ) 


treated briefly by Lin (ref. 24) and a correspond- 


ing broken-line approximation. Lin's result yields ocgZ = Ti'\pS for the 
sine profile. Equation (37) in the following analysis was used to ob- 
tain the dimensions of a profile having the same Osl. The following 

two profiles are given for flow between parallel walls, v* vanishing at 
the walls. 
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Because the widths of the profiles composed of broken-line segments 
etre so nearly equal to those ^(having identical neutral wave numbers) 
whose velocity gradients are continuous^ it is believed that the broken- 
line profile as used here has not involved the introduction of extraneous 
solutions . 


Interpreting the Results of the Stability Analysis 

The stability analysis gives an amplification rate oCj_ and a phase 
velocity Cp for disturbances as functions of the wave number a and 
the fraction of duct occupied by flame for any one station along 

the a^is of the duct. From this information a stability map is con- 
structed. The stability map consists of a plot of contours of constant 
acj^ on an oh]_ - hg/h]^ grid. This is a generalized map and contains 

no assumption regarding the rate of change of flame width with distance. 
Lines of constant frequency are shown on the stability map calculated 
from the phase velocity Cp. From the map and an assumption of the na- 
ture of flame -spreading rates ^ the total amplification of a disturbance 
of given frequency is calculated. 

Some aspects of this information are reduced to measurable quantities 
that can subsequently be compared with experimental results. These quan- 
tities include 

(1) Qualitative dependence of turbulent flame-spreading rate on 
disturbance frequency 

(2) Relation between critical wave number and flame width 

(3) Qualitative behavior of amplified disturbances as affected by 
the steady-state condition (velocity^ fuel-air ratio) 

The section Effect of Flameholder deals with the effect of a flame - 
holder of finite size on the steady-state profile and^ hence^ on the 
stability analysis. 


STABILITY ANALYSIS 
Steady-State Solution 

It might be well first to tabulate the equations derived by Tsien 
(ref. 22) relating ^2^^! fraction burned F. 

These relations are' based on the assmption that the velocity profile can 
be approximated by straight-line segments. 
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For the Incompressible case: 
Apex velocity: 



Hot-gas velocity distribution: 



Flame width: 



Fraction bTxrned: 



The hot-gas distribution (eq. (2)) and the relation between hg/h]^ 
and U]_/Uq (eq. (3)) are used in interpreting the results of the stabil- 
ity analysis. In the analysis the hot-gas velocity distribution (eq. 
(2)) is represented by 

U = Ui + (hg - y)p (y > 0) 

U = Ui + (hg + y)3 (y < 0) 

vhere 

p = (U* - UiVhg 

Nothing is said about the distribution of hg/h-j^ with distance x. 
This must come from further assumptions regarding flame slope^ flame 



speed, or rate of change of F with x. For the present it is suffi- 
cient to say that the rate of change of h 2 /hj^ with x is smallj that 
is, 

ox hi 

where X is the wavelength of a- disturbance of interest. 

The stability of the profile assumed by Tsien, which closely approx-' 
imates that calculated by Scurlock, is treated as a parallel-flow 
problem. 


First-Order Disturbance Problem 

Two excellent summary papers eire available on the stability of 
parallel flows, one by Lin (ref. 24) and one by Pillow (ref. 25). No 
attempt to review the field is made herein. The method employed here 
dates back to some of Rayleigh's early work (ref. 17) with very slight 
modification. The development of the disturbance equation that follows 
is from Goldstein (ref. 3l) . 

The flow is assumed parallel, inviscid, incompressible, plane, and 
non-heat-conducting. Gravity is neglected. The equations of motion are 

^p (hu , Su , Su\ /d\ 

S = -P (st + “ SJ 5?) 


The fluid 
with y ) , 


Sp / 

¥ = -"V 


Sv Sv , S 
§t + ^ 


Sv\ 


( 6 ) 


is assumed incompressible (but the density is allowed to vary 
so that 


Sp 

St 


Sp 

5x 


Sp 

5y 


( 7 ) 


and the equation of continuity is 
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Let 


u = U(y) + u»(x,y,t) 

V = v*(x,y,t) 
p = P + p* (x,y,t) 

P = p(y) + P* (x,y,t) j 


( 9 ) 


where primed quantities are of first order of smallness. Assume u’, 

V*, p*, and p* to vary with x and t as where a is 

real and c, in general, complex. 


By substituting equation (9) in equations (5) to (8) and retaining 
terms of the first order, equations (5) to (8) become 


ia(U - c)u* + ^ v'J = -io-P’ 

(10) 

p [ia(U - c)v*] 

(11) 

ia(U - c)p» + v*^ = 0 
dy 

(12) 

+ iou‘ = 0 

(13) 


Eliminating p', p*, and u* from equations (lO) to (13) yields an equa- 
tion for v‘ 


(U - c) 




1 ^ 
p dy 



c) 


Sv* 

5T“ 



= 0 


(14) 


It is apparent that, though equations (?) and (l2) are useful in 
defining the type of incompressible flow field assumed, they are not re- 
quired to obtain the disturbance equation (14) unless it becomes neces- 
sary to specify p' . 


Defining 


= -ia<p(y)ei°'(^-c‘t) 


V* 
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equation (14) becomes 
(U - c) 


- <p ^ i ^ f(U - c) - (P ^1 = 0 
\hy^ J dy2 P dy ^^57^ dyj 


(15) 


In regions where p and du/dy are constant, equation (15) reduces to 

(16) 




—-i - a2<p = 0 
oy^ 


with the solution 


(p = A'Le'^y + Age"^'^ (17) 

The disturbance velocities given by definition and equation (13) 

are 


u* = (A]^oe°'^ - A 2 <xe“°'^)e^°'(^“^^^ 

V* = -ia(A-Le°'y + Age"°'y)ed°'(^“^^) 


where the real part is taken in the end. 


(18) 


Boundary Conditions 

The geometry of the flow field is shown in figure 1 . The quantity 
3 is given by 3 = (U* - U]^)/h 2 . 

Equation (16) and the solutions (eqs. (18)) apply to each of the 
regions 1 to 4 considered separately. The establishment of the boundary- 
value problem consists in assigning suitable boundary conditions at the 
walls and at the boundaries of regions 1 to 4. These boundary conditions 
are; 


(1) V* vanishes at y = ± h]_ 

(2) V* is continuous at interfaces 

(3) p’ is continuous at interfaces (or equivalently, the complete 
equation (14 ) Integrated across each interface is satisfied) 
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Let the equation of the interface at y = hj^ be 


y = hi + h(x,t) 


(19) 


(in the following derivation the subscripts i and j represent condi- 
tions evaluated in a region i or j . ) Then, to the first order of 
approximation, at the interface. 


V' 


Sh 

= 5t 


+ u 


Sh 

5^ 


( 20 ) 


where U is the mean velocity at the interface. In order for the pres- 
sure to he continuous at the interface^ it is necessary that 


<^pj _ ^ 

dxi “ dx^^ 


\ 


where x^ is along the Interface, or equivalently in the x,y-system 


^Pi ^Pi hh 


^P 

57" 


+ 


^h 
dy 5x 


where at the interface 


dy _ dh 
3x “ 5x 

From the equation of motion, 

dp / . du . du\ 

dx = -P Vdt ^ 57/ 

dp . dv , dv\ 

57“ + ""d3^+ '"57/ 

where 

u = U + u» 

V = V* 


( 21 ) 


p = p + p’ 



NACA TN 3830 


19 


By retaining terms of 0 and first-order approximations. 


Sp ^p* 

= sr 


Sp _ dp* 

dy - dy 


■pIsF’ ^ — 


/ dv* „ dv 
■P ( 3+~ + U 3- 


or dp/dx and dp/dy are first-order terms. But, from equation (20), 
it is clear that dh/dx is of the first order so that the products 
(dp/dy) (dh/dx) are of the second order, and equation (21) becomes 

dp. dp. 

d3T = d/ (23: 


/dut du! duA fdu* du* 

Pi ^i dF" ""i 57/ [_5t^ d/ '"i (^55^ (24) 

at the interface between the regions i and j. Equation (24) results 
from a necessary condition for the pressure to be continuous at surfaces 
where p or du/dy undergo jumps. 

In order to show that this condition is equivalent to requiring 
that the integrated equation (14) be satisfied at the interface, equa- 
tion (14) is written 






(U - c) 5^ - I -p(U - c)a2v* = 0 


Integrating across the interface, equation (25) becomes 


r dvt du.] r 
\ 5F - '1 5tJ ' ['"j 


dy~ " '^j dy" 


which^ with the insertion of the continuity equation^ becomes 


r du! du.~| r 
i [_(Ui - sr + 'i = Oj [("j - 




But u' varies as so that equation (26) can be recognized 

as another way of writing equation (24). 
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The boundary conditions can he written 


v*^ = 0 at y = h^ 

(27 a) 

V* = v' 
1 2 

at y = hg 

(27b) 

^p’ 

X C, 

dx ** 5x 

at y = hg 

(27c) 

v’ = V* 

at y = 0 

(27d) 


sr = 53T 

at y = 0 

(27e) 

^3 = 

at y = -hg 

(27f) 

^p 4 ^Pi 

^ = 'sr 

at y = -h 2 

(27g) 

v| = 0 at y = -h^ 

(27h) 


where the subscripts 1 to 4 denote flow regions. 


Solution of the Boundary Value Problem 


The solutions of equation (l 6 ) for regions 1 to 4 (fig, l) are 
written as follows: 


Region 1 = A-i^e^'^ + Age"®'^ ^ 

Region 2 ^2 ~ 

> 

Region 3 T3 = C^^e^-y + 


(28) 


Region 4 *>*4 = D-|^e“'y + Dge"^'^ J 

Inserting equations (28) into equations (27) results in eight equa- 
tions in the eight unknowns A^, Ag, . . .Eg. The requirement that 

Ai, Ag, . . . Dg have nontrivial solutions is 



NACA TN 3830 


21 




oil 
0 

0 0 0 

0 0 .0 

0 0 0 


0 0 0 

0 0 0 

0 0 0 

-1 -1 -1 

0 



0 0 ri 


. I - p-c) - I -fu*.c) . g 



Where and ^2 = gy assigning = 0 (this 

merely alters the real part of c) and operating on rows 1^ 2 , 3^ 6, 7, 
and 8, equation (29) can be reduced to the determinant 



n _ 
m 



22 


NACA TN 3830 


On expanding and simplifying, equation (30) "becomes 


p/n Pl^ ' 

c - £-( + 

aim p„ 


where 


-1 


2 /n Pi _ n Pi n Pi \ hp2 p2 

li+1 - C hpP — +hP p)+ -i— |i - i-;;r 

_ \m P2 J m a ^ m P2 a^_ 


= 0 
(31) 


P = 


1 + ri 
^-ri 


One root can "be seen to "be wholly real 


c 


P n 


Pi 


a I m p„ n 


(32) 


This can he identified as the root belonging to a symmetric disturbance 
by assigning = 0 at . y = 0 in equation (27). In this case 

equation (30) reduces to 


c 


P + E ^ 

a m Pg 


1 



(33) 


which yields as before 


c 



-1 



( 34 ) 
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can similarly be identified as roots belonging to an antisymmetric dis- 
turbance. One way of doing this is to set at y = hg equal to 

at y = -hg or (from eq. (28)) 

+ Vi - Vz - =2 - ° (5=) 


Using equation (35) in place of equations (27f), (27g), and (27h), the 
determinantal equation becomes 



(36) 


with roots of equation (34) for c as before. That the root for the 

- / n- \ “1 


symmetric disturbance c 


Pg / 


is wholly real implies neutral 


stability. For, the symmetric disturbance^ then^ the inviscid analysis 
does not yield information on wave numbers where amplification occurs^ 
but merely suggests that the profile is relatively stable to symmetric 
disturbances and neutrally stable in the. limit of infinite Reynolds 
number. 
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The complex root having the positive imaginary part corresponds to 
an amplified disturbance. The denominator of the right side of equation 


(34), 



when written in full, 


2 



-2ohT 

+ e 

-2ahn 

- e 


-2ohp 
1 + e ^ 

-2ohn 
1 - e ^ 


+ 


is positive for that is, for the flame zone contained within 

the duct. Therefore, the amplified root of equation (34) is 



(38) 


Because this is ah inviscid analysis, the damped root of equation (34) 
can have no meaning (ref. 24) and is subsequently ignored. 


RESULTS 


CtC-st 


The disturbance amplitude grows with time as e , where acj^ 
from equation (38) is 


CLC4 = 



(39) 


Since 



AU 

h 
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The quantity under the radical in equation (39) is independent of 
U* - Ui 

local n values. This quantity is plotted in figure 2 as a func- 

“2 

tion 2 oh 2 for the following conditions: 


^2 


1 


T "^2 
V, ” hi 


= 0 


Pi ho 

-^ = 5, 7 = 0.25, 0.50, 0.75 

^2 


The first set 'of conditions corresponds to Rayleigh's solution for 
the stability of an unconfined jet. The latter set of conditions repre- 
sents density ratios that are expected in a flame in a duct for three 
flame widths. 


The values of the parameter 0Cj^(2h2/AU) in figure 2 go from 0 
at 2 oh 2 = 0 through a maximum, then to 0 again at a finite value of 
2oh2. For a given veilue of 2h2/AU, this parameter behaves as the am- 
plification rate ocj^. 


It is clear from figure 2 that, for constant density, the flow is 
more unstable the smaller h^h^ becomes. For constant h 2 /h 2 , the 

flow is more unstable the smaller P]^/P 2 becomes. 


The influence of density ratio and wall proximity on the stability 
of the subject profile is in keeping with experimental observations of 
another type of profile with a flex, namely, the wake of a cylindrical 
rod. The stabilizing effect of density ratio is illustrated in refer- 
ence 1, wherein it is noted that a flame anchored in the wake of the rod 
inhibits vortex shedding. The stabilizing effect of wall proximity is 
illustrated in reference 32 wherein it is noted that wall proximity in- 
creases the minimum Reynolds number for the occurrence of shed vortices. 

The stabilizing influence of the walls is further illustrated in 
figure 5 . Figures 3 fa) and fb) show the values of 2odi2 for neutral 

stability and maximum instability plotted against h 2 /h 2 ^. Figure 3(c) 
shows the amplification parameter ocj^( 2 h 2 /AU) for maximum amplification 
rate plotted against h 2 /h]^. Both the amplification rate and the range 
of the parameter 2 oh 2 diminish as the flame fills the duct. 
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Construction of Stability Map 


In considering the effect of disturbance amplification in the flame 
zone, local-growth-rate contours as a function of wave number and flame - 
width ratio ^ 2 /^! Present nearly all the important information. Such 

a plot is called the stability map of a flame in a duct. The map is con- 
structed from the data of figures 2 and 5 and from the steady-state re- 
lation of U*/Uo functions of hg/h]^ (eq,s. (l) and (3)) 

Ui /hgX u* A2\ 

in the following manner: Given: I hw^ U” I hw^ 

2hp / hp\ Pi 

oci I 2oh2, 1 for a given — . The quantity desired is 

2hi / hgX 

OC 4 ^ — l2ohn, ^1. This quantity is obtained through suitable cross 
^0 \ “ 1 / 

plots of fij^es 2 and 3 and equations (l) and (3) by setting 
2ohi = Zah^f (h2/hi) : 


yl - lii 

2hi 2h2 Uq Uq 

Uq “ '“i U* - Ui "EiThi" 


(40) 


The resulting stability map is shown in figure 4 . Contours of con- 

2hi 

stant dimensionless amplification rates oci — are shown as functions 

^0 

of the dimensionless wave number 2ohi and dimensionless flame width 
h 2 /h^ for Pi/p 2 = 7 (fig. 4(a)) and Pi/p 2 = 5 (fig. 4(b)) . These con- 
tours are in a region bounded by 2ohi = *, 0^ ^ 2/^1 ~ line 

of 2oh-j^ for neutral stability, that is, the critical 2oh-|^ as a func- 
tion of h 2 /h-L. 


Phase Velocity 

A disturbance that distorts the flow and therefore distorts the flame 
front has a wavelength X = 2jt/a. Such a disturbance if introduced at 
the flameholder travels along the flame front at the phase velocity Cj.. 

This can be seen from the variation of the disturbance with x distance 
and time; that is. 


ia(x-ct) otc,-t ia(x-c_t) 
v'~e =e-‘-e ^ 
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In order to follow the amplification rate of a given disturbance, it is 
necessary to follow a disturbance of constant frequency. 


The relation between frequency and wavelength is 


f 


Cr CLCj. 


(41) 


The quantity 


^r - ^1 
U* - Ui 


from equation (38) is given by 


Cr - Ul 

U* - Ui 


ahc 


n ^1 

^2°- ^ - 1) 


for 


[“ 2 “ * S ^ (>“*2 ■ ^ ^ ‘ ° 


and by 


n Pi. 


Cr-Ui 1 h2a+-— (4oh2-l)+ 

U * - Ui °^2 


(42) 






(43) 


for 


^ (^^°^2 (°^ 2 P - 1 ) 2 : 0 


m p, 


m p.. 


A plot of this quantity against oh 2 is given for P]^/p 2 = 7 at 
^ 2/^1 = 0*25, 0.50, and 0.75 in figure 5(a) . Also shown in figure 5(a) 
for comparison is (c^. - U 3 ^)/(U* - U^) for Pi/p 2 = 5, ^ 2/^1 = 0*25 and 

Pj^/ Pg — 1, ^2^^1 ~ ® * 
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Note that 


lim 

ohg"^ 


^r-U l 

U*-U 


From the steady-state solutions, 

lim (U* - Ui) = 0 

hg/hi-^O 


Therefore, at the flameholder Cj. = = Uq* Then it is convenient to 

relate the local wave nvuniber to the wave number at the flameholder by 


a = ttr 


cr 


(44) 


where Oq is the wave number at the flameholder. It is related to the 
frequency of the disturbance f and Uq by Oq = f2n/UQ. By means of 

equations (42), (43), (44), (l), and (3), values of a/o^, ^ were cal- 


culated and they are shown in figure 4 as lines of constant frequency. 
These lines enable one to follow the local amplification rate of a dis- 
tiirbance of given frequency from the flameholder to some station down- 
stream and are so employed in a later section. 


An additional calculation for the phase velocity was made for later 
comparison with experimental data. In figure 5(b) is shown a plot of 
Cj,/u^ against h^/h-j^ for “ 7 and 2cxQh^ = 8, 12, 16. These data 

enable one to compare directly measured values of phase velocity with 
those predicted from the stability analysis. 


Disturbance Growth 

Figure 6 presents a qualitative picture of the nature of the dis- 
turbances propagating at phase velocity. 

The streamline distributions are given for a neutral disturbance in 
a reference plane moving at the phase velocity Cj.; in this frame of ref- 
erence, the streamlines are no longer a function of time. There is mo- 
tion relative to the "cat's eyes" of the disturbance centers. In the hot 
flame core, the gas moves faster than c^.; in the colder gas outside the 

flame, the gas moves more slowly than Cj.. 
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To picture what happens to an amplified disturbance, imagine the 
y-amplitude of the wrinkles in figure 6 to be growing from right to left. 
In addition the flame width is increasing gradually^ c is also increas- 
ing gradually, and, therefore, the x- and y-spacing of the "cat's eyes" 
is increasing in the direction of flow until a neutral situation is 
reached. 


At each x-station the disturbance has a local growth rate acj^ that 

acts over a short time At. Then (since x is related to time by 

dx = Cpdt) at any time t after the distvirbance originated at the flame- 

holder, the amplitude is 


v» 


Pt 


cxcj^dt 


(45) 


where is the amplitude of the disturhance at the flameholder on the 

2 

centerline of the profile. 


Up to this pointy the independent vaxiahles have been the dimension- 
less wave number and flame width h2/h^. In order to evaluate 

equation (45)^ it is necessary to use the relation 


dx 

dt 


^r 


(46) 


and some assumption concerning the variation of h^/h^^ with x. Two 
simple assumptions are used that bracket the variation of ^ 2/^1 

X found by experiment and facilitate the integration of equation (45). 
These assumptions are (l) that the flame is linear: 





(47) 


and (2) that the change in fraction burned with distance is linear: 

= Kg (48) 

d ^ 

^1 


^The linking of the parallel flow solutions by conserving v* along 
the centerline is the simplest but not necessarily the best of several 
available choices . 
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A qualitative comparison of the way these assumptions describe a flame 
svirface is given in the following sketch: 



Inserting equations (46) and (47) 


in equation (45) results in 


v» 





(49) 


The term in parentheses is called the gain g. It is independent of the 
value of the propagation rate 


A similar derivation can he made for the alternate form of heat 
release^ that is^ for 



vhich yields 


V' = 



(XC 


2hi Uq 


,1/K2 


1 Ur 


— dF 


(50) 


2 - 


where the term in parentheses g* 


is independent of the value of K' 
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The values of g and g* can be obtained by graphical integration. 
They give relative values of amplification of a disturbance of given 
frequency when integrated along a constant-frequency line in figure 4 . 

In figure 7 are given the results of such an integration carried 

from hg/h-^ = 0 to hg/hj^ = 1 for a range of 2aQh-j^ and for the follow- 
ing conditions: In figure 7 (a) for Pp/P2 = 5, 7, g is given as a func- 
tion of 2aQh^. In figure 7(b) for Pp/ Pg = 7, g* is given as a function 

of 2oQh-]_. It is apparent from figure 7 that, for the conditions examined, 
there is a most -amplified frequency corresponding to a value of 2oQh2^*6. 
Recall f = 2aQh-j_UQ/2n2h2_ or the most-amplified frequency is approxi- 
mately f ® lJQ/2h^. This frequency does not depend upon a rate of flame 

spreading. It was obtained by assuming that a mode of flame spreading is 
specified. Approximately the same frequency obtains for the two modes 
of flame spreading specified. 

In figure 7(a) not only is the most -amplified frequency at P]^/Pg = 5 

very near that, for P]^/Pg = 7, but the numerical values of the gains are 

very nearly the same as well. This latter occurrence is’ due to two oppos- 
ing influences: For a given U* - U-j_, the profile having the lower value 

of P]^/pg has a higher amplification- rate, and a disturbance of given 
frequency is amplified over a wider range of hg/h^^. On the other hand, 
at any value of hg/h^, U* - U-]_ is smaller for the smaller p|j^/pg. 

These opposing influences yield the results shown in figure 7(a) . 

The growth of the disturbances for fractional values of hg/h;j_ (or 

F) was calculated by taking the integral in equation (49) or (50 ) to 
fractional values of hg/h]_ (or F). The results are shown in figures 

8(a) and (b) respectively, in which the gains g and g‘ are plotted 
against hg/h^_ or F for values of 2aQh-]_ from 2 to 18. It can be 

seen that, for some values of hg/h-^ < 5, there are frequencies higher 
than those corresponding to 2oQh]^ ~ 6 that locally exhibit a maximum 
gain. The fact that for values of hg/h-|_ <0.5 there is a local ampli- 
fication that depends on frequency is used later in comparing analysis 
with experiment. It can be seen that as hg/h^^ decreases, the frequency 

that is most amplified increases. 


Disturbance Velocity Distribution in y-Direction 

The way the disturbance velocity v‘ varies along the axis of the 
duct is shown in figure 8 . In reducing the results of this analysis to 
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measurable quantities, there is need for a qualitative picture of the way 
the disturbance velocity v* varies across the duct. Such a picture can 
most easily be obtained by considering a neutrally stable disturbance. 

From equations (18) and (28), 



h 

Be 


+ rf 


^y n 

^ + 1 

^2 


0 ^ y ^ h 2 


(51) 


where 


2,y 




^y " 


disturbance velpcity at y in region 2 
disturbance velocity at centerline 
e-ay 


In the region h^ < y < h^, v£ y/v2,h2 obtained from equations 

(18), (27), and (28) in the form 


V- 

V. 


I ^2 . „2 Y- 

l,y ^1 ^ ^y ^ 


ih2 -rl + rl 


(52) 


The ratio B]_/B 2 can be obtained from the equation represented by 
the first row in the determinantal equation (30) 




+ 1-c 


a • m pg ° 


or 


Be 


_r2 

^2 


/E: ^ _ I'j _ J- 

r P2 / ^ 

L\“ P2 / 


( 53 ) 
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Recall that p/<xc in equation (30) is the same as 


1 u*-ui 

<xh2 C - 


> which 


for real values of c is given in equations (42) and (43), 


A series of calculations was made to obtain the representative be- 
havior of across the duct, and the results are shown in figure 

^2 0 

^ for p-^/ p2 = 7 . The critical h2/h-|^ is taken as the value of 

where for a given frequency the imaginary part of c becomes 0. Figure 
9(a) shows the value of v^ /v^ at the critical hg/h^ plotted against 

2 

initial wave number 2aQh-|^. It can be seen from this figure that for 
2oohi > 8 (i.e., for (h2/hi) critical < 0*^6) the value of vfjg/v^ • 0.2. 
The ratio falls off as 2oQh2^ < 8 (i.e., for (112/h]^) critical ^ 0.46) as 
the flame approaches the wall. 

In figure 9(b) the value of the ratio is given for a fixed 

2 

frequency corresponding to 2oQh-|^ = 12 following a disturbance along the 
duct from the flameholder (h2/hj^ = O) to the station where the flame 
touches the wall h2/h-L. The dotted portion is interpolated between 

\/\ = 


0 where v* 


= v; 


and hg/h^ = (hg/h^)^ 


where the cal- 


culation starts, 


hg ’ 0 “2' **1 *‘l ' critical 

It is obvious from this curve that^ if one measures 


h2 


instead of v^, the apparent growth rate is not as rapid for 

*^2/^1 (^2^1 ) critical because of the velocity dis- 

tribution. Also, for h2/h2^ > (h2/hi)critical ’there is an apparent 
decay in v^^ owing to this distribution, despite the fact that the 

value of v' is neutrally stable. This latter point is employed in com- 


paring the analytical with the experimental results, 


The significance of figure 9 is summarized as follows: 

(1) Although the inviscid treatment eliminates consideration of 
damping, it is clear that, at the flame front, the velocity amplitude 
diminishes at h2/h-L > (^2/^1) critical h^cause of the velocity 
distribution. 


(2) In some instances, even though the flow is unstable and the 
disturbance velocity amplitude is increasing along the center of the 
flame, the distiirbance velocity amplitude at the flame front may not nec- 
essarily be increasing. 
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(3) The disturbance that is most -amplified at the center of the 
profile (i.e., 2oQh-j_ = 6) is not necessarily the most-amplified at 

the flame front, owing to the velocity distribution; the frequency of the 
latter is somewhat greater than the frequency of the former. 


Relation of Theoretical Results to Measurable Quantities 

In the experimental phase, current instriomentation techniques im- 
pose limitations on what can be measured. For example, the velocity v^ 

cannot be measured directly by hot-wire techniques because the flame, tem- 
perature is too high (»4000^ R). The purpose of this section is to col- 
lect certain results of the foregoing analysis that lend themselves to 
direct measurement. 


Flame slope . - In equations (49) and (50) the growth of a disturb- 
ance v*/v^ is found to be proportional to a calculable quantity, gain. 


hi 


raised to a power 1 /k. The quantity K is alternately K]^ = or 


^1 


Ko = 


dF 

^1 


depending on the mode of flame propagation assumed. 


For both modes, K is smaller the more gradual the rate of flame 
spreading. To a first approximation, then, the amplification of a dis- 
turbance increases as the ratio of propagation velocity to flow velocity 
decreases. This trend can be checked experimentally. 

Phase velocity . - One quantity that can be measured immediately is 
the phase velocity c^. A plot of c^/Uj_ for three frequencies is given 

figure 5(b) . Thus, for a given flow setting both the absolute magni- 
tudes of Cp/Uj^ can be compared and the trend with frequency and with 

h^/hf variation can be examined. 

Flame-front displacement . - The transverse velocity v* acting at 
the flame surface y = h2 causes a wavy wrinkle to appear in the surface. 

If a quantity h is defined as the distance from crest to valley of these 
wrinkles, this quantity equals twice the y amplitude of the streamlines 
of the flow at the flame front. If a neutral case is considered and the 
distribution of v* in the flame zone is approximated by 


+ (»4 - ^.)(i - 
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the equations of the streamlines are given by 


— *» 1 - 6 cos ax i 

yc 




2Y + 2K cos <xx 


where 

coordinate where U = c 


^ - '"c 


(U - c)ay 


¥* dimensionless stream function. 


s - 


(U* - o)y<, “'o 


Udy 


K 


(U* - c)ay^, 

v^ amplitude of disturbance velocity at centerline 

v^ amplitude of disturbance velocity at y^. 

At hg ^ hg critical^ varies between 0.96 and 0.99. Because 2K 

be taken arbitrarily small, for 2K/(1 - 2Y^ <<1^ 


can 


yc 

or, for y > Yq, 


— ?= 1 - 6 cos OX ± 




1 + 


K cos ox 


[_ (1 - 21 *) 


1 


-^ = 1 + Vl - 2Y* + 
yc ^ 


K 


Vl - 2f^ 


6 \ COS ax 


1 + 0.2 + 


a-yr 


v^/Uq '^hho - '^c/Uq 


ir - 


l_ «0 


u - u. 


0 


Un -J 


«» 1.2 + 


o-y, 


^6/Uo 

u*- 

C Un 
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But v^/vq is a weaJt function of a and hg/h-i^ and is approximately 

0 . 2 . 

Employing the identity 

U - “l/ 

and a value of (c - Up)/(U* - U]^) from figure 5(a) 



£ 

Ur 




(0.85) 


where the quantity U*/Uq - U]^/Uq is a function of Pi/p2^ only. 

Without going into further detail, then, some predictions can be made 
from the theory that may be checked by experiment. First, h may be ex- 
pected to behave qualitatively as shown in figure 10. 


The shape of the curve in figure 10 can be used to locate the ap- 
proximate X position of There is a point of inflec- 
tion where ''^h 2 /hi maxinMm and a knee as v 

to grow. The knee position falls between the h 2 /h^ for maximum ampli- 
fication and (h2/hi)critical* Thus^ the data that can be sought are the 
absolute values of h 2 /h^ where the knee occurs for a given frequency^ 
and the relative positions of the knee for different frequencies. 


h2/hl 


ceases 


Flame propagation . - As stated in the HJTR0DUCTI0N , one of the rea- 
sons for pursuing the stability analysis is to find a disturbance fre- 
quency that permits the greatest increase in flame spreading with the 
smallest expenditure of energy. In this paragraph a model of flame prop- 
agation is constructed that permits extrapolating the behavior of v' 
from the stability analysis to a measurable behavior of flame propagation. 


The flame front is assumed to be a zone in which the density goes 
from to pg in a distance that is short compared to the disturbance 

wavelength. This front is propagating into the gas of density p^ at a 

rate that is, the undisturbed turbulent flame speed. Owing to the 

normal propagation of the front, an initially wavy surface becomes a 
cusped surface. Assuming such a surface composed of cylindrical arcs 
(following sketch) 
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permits a calculation of the flame-front stretching due to a displace- 
ment h. If the undisturbed area Ay = X and the distvirbed area 

Ay are approximately equal to arc length. 



(54) 


An equilibrium h is reached when 



(55) 


Then Ay/Ay behaves roughly as fl + ^ ^ — ) for small amplitude. One 

\ ^ ^,U/ 

sees that equation (54) contains a nonlinear term 4h^/X^. The develop- 
ment so far has been a linear one. The linear development would be ex- 
pected to furnish only a first approximation for amplitudes large enough 
so that 4h^/X^ is no longer negligible compared to unity. 


The sequence of events that leads to a prediction of Ay/Ay is: 


(1) The flame is initially wrinkled by the velocity v* acting at 
the flame surface. 


(2) As the wrinkles grow, the flame shape shifts from a sine wave to 
a cusped surface, owing to normal propagation. 

(3) After a certain time an equilibrium is reached between the dis- 
turbance velocity v* and the undisturbed propagation velocity y. 

The ratio of the disturbed surface to the undisturbed surface when 
equilibrium is reached is, for small disturbances. 


^ = (l + 1. V' ^ ^ ^T,D 

% \ 2 Srp^yy S|j,^y 


( 56 ) 
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Before this equilihrlum value of Ap/Ay is reached, the flame surface 
is given by (again for the neutral case) 


yc 


0.2 vA 

« 1.2 + 77 = COS ax 

AU ay^. 


But 2h/X of equation (54) is given by 


2h 

X 






yc2a 

Jt 


AU jt 


cos 


ox 


Thus, in the transient case as in the equilibrium case (eq. (55)) Ap/A^ 
behaves as v* . The quantity Ap/Ap when measured locally should behave 

like g in figure 8 . if the effect of distribution is neglected. When 
the latter effect is taken into account, v’ for fixed Soph-L measiared 

at h 2 achieves a maximum value somewhere between hg/h^^ for maximum 
amplification and neutral stability. 

Thus, at any station, hg/hj^, there should be a frequency of disturb- 
ance that would cause the greatest increase in flame speed and, at this 
station, the increase in flame speed should be greater for greater values 
of v’/s^^p. 

Figure 11 gives a plot of 2ophp against hg/hp for maximum ampli- 
fication and for neutral stability. For any value of hg/hp, the fre- 
quency giving the greatest ''^'/®T,U flame surface) has a value 

of 2aphp that falls between these curves. The lag in time required for 

the flame to become wrinkled initially will move the curve for maximum 
Ap/Ap to slightly higher 2ophp values. The curves given should be 

considered the lower limits of 2aphp for maximum Ap/Ap, and they serve 
as a basis for comparison with the experimental data. 


Effect of Flameholder 

When a flame is anchored in a duct, some form of flameholder must be 
used. The most common form of flameholder studied in recent years is the 
cylinder with axis normal to the stream. The cross section of the cylinder 
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has assumed many different shapes. There are properties of the flow in 
the waJce common to all the shapes studied. In this section some of these 
properties are used to deduce the effect of a flameholder on the flow 
field far from the flameholder and, hence, on the stability of the flow 
field far from the flameholder. 

Steady-state velocity distribution . - The properties of the flow 
field in the imm ediate wake of the flameholder as described in references 
1 and 33 are shown in figure 12 . There. is an attached eddy pair in the 
wake of the flameholder. The gas in this zone is hot, acquires energy 
from the downstream hot and delivers energy to the outside cold gas 

by diffusion. 

At some position on the surface of the recirculation zone, suffi- 
cient energy is acquired by the cold gas for ignition to occur and a 
flame to propagate into the cold-gas stream as shown in figure 12 . The 
pressure measured at the walls has an axial distribution shown qualita- 
tively in figure 12 (refs. 1 and 33). 

If the static pressure is assumed constant across the duct at any 
ayj station, the effect of the first positive, then vanishing, then 
negative pressure gradient upon gases of unlike densities permits an 
estimation of the velocity profile downstream of the flameholder. 

In figure 13 are sketched the velocity profiles at stations a, b, 
and c of figure 12 . At station a the pressvire gradient in the axial 
direction has been positive for a small distance and is still positive 
at a. This fixes the sense of the velocity difference between cold and 
hot gases on either side of the flame and the sense of the velocity 
gradient of the hot gas at the flame front. At station b the gas has 
had a positive pressure gradient history, again fixing the sense of the 
velocity difference, and at b the pressure gradient vanishes, causing 
the velocity gradient at the flame front to vanish. At station c the 
gas has a negative pressure gradient history from b sufficient to cause 
the sense of velocity difference and gradient at the flame front as shown 
in figure 13(c) . The gas particle at the position of the velocity maxi- 
mum passed through the flame at station b. If the prpfile is approxi- 
mated by straight-line segments with a corner at this velocity maximum 
as shown in figure 14 , some quantitative relation based on the continuity 
equation and experimental observations can be obtained. 

From continuity. 


1 * ^1 

i b«U* = - 


( 57 ) 
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and 


U* + Ui Pn 

b = Sflf - 


(58) 


where 


p /p„ density ratio of cold to hot gas 

cold-gas velocity 

U* maximum hot-gas velocity 

b,b' defined by fig. 14 

axial length of wake to station b (fig. 12) 


If axial length from station b to station of interest downstream 

of b 

mean flame speed over 
Sf mean flame speed over If 


Combining equations (47) and (58) gives 

IV + Pj Vw 

2 b 2U* Sflf 


(59) 


1 b' 

In order to obtain an order of magnitude of 2 ~ interest in the 

stability analysis, the terms on the right side of equation (59) are 
given the following values: (U* +U]^)/U* = 1.5 (from ref. 1 for F = 0.25, 


Pj^/Pg = 6), where F 


is fraction burned; 


Sf Sqi/Sl 


1 

2 . 5 ^ 


where 


Sgn/Si 


is ratio of turbulent to laminar flame speed; » 4D (from ref. l), 
where D is flameholder width. Defining I as + If , 


I 



1 


4 



1 
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equation (59) becomes 



At a station Z/D = 8, from equation (60), 




Farther downstream, at z/d = 24 for example. 


(60) 



The foregoing stability analysis is for a profile in which b* = 0. 

From this section it can be seen that, near the flameholder 
(Z/D =8), the hot-gas velocity distribution differs in shape from the 
profile examined neglecting the flameholder. At this station 60 
percent of the flame-zone width has a flat distribution. Remote from 
the flameholder (z/D = 24), 6 percent of the hot-gas width has a flat 
velocity distribution. 

In the next section the magnitude and direction of the error in- 
curred in neglecting this flattening is estimated. 

Comparison of the stability of a profile and a -7"^ profile. 

2 a comparison is made of amplification rates as a function of 
2 oh 2 for several values of P]_/p 2 seen that the 

effect of increasing P]_/p 2 stabilizing; that is, the 

magnitude of the amplification rates and the range- of 2 oh 2 over which 

positive values exist are reduced. The plots of amplification rate 
against Sohg are affine. It appears reasonable that, in examining the 

effect of truncating the profile, an examination at a single value of 
P 1 /P 2 and h^/h^ would yield representative results. 


The conditions considered here are the ones formulated by Rayleigh 
in equations (24) and (25) (p. 387, ref. 17) and correspond to the case 
where P]^/ P 2 = ^ ^2^^! ” ^ h^ = «) . Calculations were made 

1 b’ 

in which ^ ^ was taken as a parameter in such a way that the results 
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yield values of 00 ^ — and 2ab at maixinnim instability and a critical 
value of 2ab where ocj^ vanishes. The results of these calculations 
are shown in figure 15 . 

When b* = 0^ the results correspond to those of figure 2; the flow 
is unstable to antisynnnetric disturbances, neutrally stable to symmetric 
disturbances. 

1 b' 

As 2 ~ increases, the amplification rate and the critical values 

of 2ab decrease for the antisymmetric disturbances. The flow gradually 

1 b' 

becomes unstable to symmetric disturbances as 2 ~ increases. At 

^ ^ >2.0, the flow is as unstable to symmetric as to antisymmetric dis- 
turbances and possesses the same critical values of 2ab. The values of 

1 b * ■ 

— — of interest here are of the order 0.6 to 0.06. The flame half 

2 b 

width h .2 is i b* + b. In the following table the values of the im- 
portant parameters clc^ ( 2 ah 2 )^g^ and (^ohg) critical 

1 b* 

pared for 2 “ = 0.06, 0.6 for antisymmetric disturbances: 


1 b* 

0 

0.06 


2 b 

0.6 


0.492 

0.508 

0.688 

(20-^2 )max 

2.44 

2.49 

3.28 

(2c‘^ 2) critical 

3.69 

3.71 

4.57 


For a given 2h2/AU, the amplification rate is increased as the hot-gas 
velocity is progressively truncated. 

The wave-niimber range over which the flow is unstable is widened. 
Both of these effects bespeak a greater instability of flow.' On the 
other hand it can be easily demonstrated that, for a given AU, h 2 /h-j_ is 

greater when the profile is truncated so that the destabilizing effect 
of the flameholder is partially but not completely attenuated. 
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1 b* 

The difference between the parameters calculated at ^ ~ ® 

0.6 shows that near the flameholder the amplification rates may be of the 
order of 20 percent higher for some wave numbers than were those previ- 
ously found. Because the flame zone must start with a finite width de- 
pendent upon the flameholder width (say, D), those frequencies having 
such wave numbers that 2oDK exceeds are suppressed. 

Here, K is some proportionality factor less than unity. 

Certainly, neglect of the effect of a finite-size flameholder facil- 
itates generalizing the results of the stability analysis. The investi- 
gations of this section indicate that for flameholders small as compared 
to duct width, there is a large part of the flame zone in which such . 
neglect appears warranted. 


SUMMARY OF THEORETICAL RESULTS 

1. The flow field of a flame in a duct is unstable to antisymmetric 
disturbances having wavelengths of the order of 2.5 times the local flame 
width and greater. As h 2 /h]_ approaches unity, this critical wavelength 
approaches 0. 

2. The amplification rate depends directly upon the magnitude of the 
velocity difference (apex velocity minus cold-gas velocity) and inversely 
upon the flame width. For fixed values of these quantities, the maximum 
amplification rate decreases with increasing values of P-,/pp arid h /h , 

going to zero for hg/h-]^ = 1. 

3. When the history of a disturbance is followed through the course 
of an entire flame, it is found that a wave number for maximum disturb- 
ance growth obtains. This maximum growth is for such a wave number at 
the flameholder that 2(XQh]^ = 6 which corresponds to a frequency 

f = UQ/2h]^. The actual growth is substantially independent of density 
ratio. This is due to two compensating factors: 

(a) The amplification rate is lower for a higher density and given 
velocity difference. 

(b) The velocity difference is greater for a higher density ratio. 

4. The disturbances achieve their terminal amplification at a value 
of h 2 /hj^ roughly inversely proportional to the frequency of the dis- 
turbance. Thus, higher frequencies achieve their amplification near the 
flameholder; the lower frequencies, remote from the flameholder. 
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5. The magnitude of the terminal amplification depends upon the in- 
verse of flame slope or roughly (Sm/Up.)"l. The amplification that can 

Uo/Sm 

be achieved for the maximum of 2(XQh^ = 6 is about 1.1 ^ For ex- 

ample, for Sqi/Uq = 1/20, amplification is 8-fold. For S^/Uq = l/lOO, 
amplification is 14,000-fold. 

6. The disturbance propagates with a phase velocity up to 60 percent 
greater than the cold-gas velocity. 

7 . The f lameholder serves to inhibit the growth of antisymmetric 
distiirbances having higher frequencies but permits symmetric disturbances 
to be amplified near the f lameholder. At 8 diameters downstream of the 
flameholder, the stability picture is substantially that obtained by . 
neglecting the flameholder. 


EXPERIMENTAL EXAMINATION OF DISTURBANCE GROWTH IN A V-FLAME 

INTRODUCTION 

The experimental research investigated the effect of imposing con- 
trolled disturbances on a V-flame anchored in a duct for two reasons: 

(1) To gain quantitative information on the mechanism of the break- 
down of the flame -generated shear region that would explain the 
occurrence of flame-amplified turbulence 

(2) To find a type of excitation that promotes flame propagation 
without causing shell failure as does combustor resonance 

The method of excitation consists in exciting a column of gas in 
the flameholder support tube. This method offers a control on the am- 
plitude when excitation is accomplished with a speaker, but could be ac- 
complished by means of an edge tone excited by the flowing stream for 
application in an engine. Only the speaker excitation is used here. The 
sound issues from a port at the flameholder and disturbs the flame at the 
flameholder. The two types of ports used operated as follows: 

(1) Sound issued from a port just upstream of one flameholder edge, 
causing the fluctuating flow about the flameholder to be quasi- 
ant i symmetric . 

(2) Sound issued from the center of the flameholder, causing the 
fluctuating flow at the flameholder to be symmetric. 

The principal disadvantage in this form of excitation is that, for. 
some frequencies, modes of burner cavity resonance are excited. An ob- 
jective in the experimental work, therefore, is to determine .whether an 
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observed flame front distortion accompanying an imposed disturbance is 
due to flame-zone instability or to some form of duct resonance. 

No effort was made to determine effects of fuel-air ratio or fuel- 
type variation. The primary variable investigated was the frequency of 
excitation and its effect on the flames studied. The only attempts at 
correlating the results are the comparison with the results of the pre- 
ceding stability analysis. 

The chief frustration encountered is the impossibility of obtaining 
a laminar V-flame at velocities where values of Sqi/Uq are sufficiently 

low for the preceding theory to be valid. The reason for this difficulty 
is the instability of the boundary layer, which separates from the flame- 
holder and undergoes transition to turbulent motion quite close to the 
point of separation. All the results are obtained with the flow already 
somewhat turbulent. 

The primary tool in measuring flame speeds is the photomultiplier 
probe, first applied in a rigorous way to the study of turbulent flame 
speeds by Clark and Bittker (ref. 34). An independent confirmation of 
the method is obtained by comparing combustion efficiency measiired by 
the photomultiplier probe and by momentum pressure drop. 


APPARATUS 

Burner Test Section 

A photograph of the burner (fig. 16) shows the important components. 
Air and propane were filtered, metered, and mixed before entering the 
two-stage calming section. In this section the mixture passed through 
six banks of 200-mesh screen and then through a 3/8-inch-radius nozzle 
plate into the l/2- by 2-inch (or alternately 1- by 2-in.) rectangular 
glass-walled test section (fig. 17) . The top and bottom walls of this 
test section were l/2- by l/2-inch (alternately l/2- by 1-in.) brass 
stock. The side walls were l/4-inch Vycor plates, resting on brackets 
and held against the l/2-inch top and bottom pieces by clamps. 

The exhaust system consisted of a 6-inch-diameter pipe connected to 
a hood and a blower to force the gases outside. A 4- to 6-inch spacing 
was maintained between the test section and the exhaust duct, and between 
the exhaust duct and the hood; no indication of resonance due to the ex- 
hauster was encountered. 

The flameholder consisted of a flat hollow paddle terminated on one 
end by a 60° gutter and on the other end by a 3/8-inch tube leading to a 
piston-type loud speaker. The three flameholders used in this program 
are shown in figure 18 . The flameholders used to supply the disturbances 
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were resonators that had resonant frequencies every 190 to 200 cycles 
when coupled to the speaker. The flameholders ranged from 36 to 38 
inches in length. The flameholder designed to produce antisymmetric 
disturbances allowed the sound to issue from a slot just upstream of one 
edge of the flameholder. The one designed to give symmetric disturbances 
allowed the sound to issue from the center. 

The flameholders pictured were approximately 0.25 to 0.30 inch thick. 
In one test a 0.50-inch gutter was tacked onto the 0.306-inch antisymmet- 
ric flameholder to examine the effect of flameholder size. 

In no case was the distribution of v’ downstream of the flame- 
holder measured. It is not known how much the experimental distribution 
differs from the measured theoretical distribution. The essential fea- 
tures in both experimental and theoretical disturbances are that anti- 
symmetric disturbances cause the flow to flap back and forth at the 
centerline of the duct^ while symmetric disturbances cause it to dilate 
about the centerline. 


Auxiliary Equipment 

Sound producing system . - The speaker was supplied a sinusoidal 
signal from a signal generator through a power amplifier. The speaker 
and flameholder could be positioned as desired. 

Flow metering . - Air and propane were metered through rotameters. 

Hot-wire equipment . - An NACA constant-temperature^ frequency- 
compensated^ hot-wire anemometer was used (ref. 35). Spectra and u* 
were measured at various locations in the duct. Measurements of u* 
were made primarily with a vacuum-tube voltmeter. Two spot checks were 
made against a root -mean-square computer ; and the correction factor was 
found to be approximately unity. Spectra were measured with a wave 
analyzer. 

Pressure measurements . - Velocity surveys were made with a 0.018- 
inch hypodermic tubing^ and alternately a 0.030-inch stainless tubing^ 
for impact tubes. Measurements were read on a water micromanometer. 
Velocity measurements were made at the plane of exhaust with room pres- 
sure as static. 

Photomultiplier probe . - The photomultiplier tube was installed in 
a probe as shown in figure 19 . Light reached the photomultiplier tube 
after passing through two vertical slots^ which were l/4- and l/l6-inch 
wide and 2 feet apart, and through filters if desired. Photomultipliers 
were supplied by a ripple-free regulated d-c source. 

Voltmeters with a-c and d-c current were used where indicated. 
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Shadowgraph system . - The optical arrangement used in taking shadow- 
graph pictures is shown in figure 20 . The light source could be used 
continuously or as a 4-microsecond flash. The camera was focused 1.5 
feet in front of the test section. The shutter and flash were manually 
synchronized. 

Stroboscope . - A stroboscope consisting of a slotted aluminum disk 
and a synchronous motor was used to view the disturbed flame in the pre- 
liminary phases of the investigation. The stroboscope can be seen in the 
lower left corner of figure 16 . 


Properties of System 

Velocity profiles upstream of flameholder . - Figure 21 gives the 
velocity profiles measured 3 and 12 inches from the nozzle; in each case 
the flameholder is positioned so that the measurement is taken upstream 
of the gutter. Measurements were also taken at the plane of the nozzle, 
and no variation greater than l/2 foot per second from 55 feet per second 
was found. These measurements are for a nominal flow setting of 50 feet 
per second. 

Flame photographs were taken and flame-speed measurements were made 
at nominal velocities of 50 and 100 feet per second in l/2- and 1-inch- 
wide test sections; the profiles of figure 21 therefore represent the 
highest percentage boundary layer of those encountered. 

Velocity profile in wake, no flame . - Figure 22 gives the velocity 
profiles in the wake of a 0.306 -inch -wide, 60° gutter at a nominal ve- 
locity of 46 feet per second. This set of profiles shows a comparison 
of the spreading of a wake with and without an antisymmetric disturbance 
at 200 cycles per second. The only amplitude measurement available for 
this disturbance is a root -mean-square u* of 5 feet per second, meas- 
ured l/8 inch above the lip of the flameholder. 

The effect of other frequencies on wake spreading was examined by 
placing total -pres sure tubes 3 inches downstream of the flameholder at 
the center of the wake and varying the frequencies. The results are in- 
terpreted in the following manner: Read the velocity while disturbances 

are being introduced. Find the distance from the wake at which this ve- 
locity is attained without the sound (say x) . An index of the lengths 
required to accomplish equal degrees of mixing without and with the dis- 
turbance is then the ratio x/3. A plot of this ratio against frequency 
is given in figure 23 . 

Because the apparent diffusion rate decreased at frequencies greater 
than 1000 cycles per second, it was thought this might imply a decrease 
in form drag for frequencies greater than 1000 cycles per second. Static 
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measurements were made at the base of the ^tter. No change in base 
pressure was observed for the higher frequency disturbances. Those lower 
frequencies which increased mixing also decreased base pressure^ that is^ 
increased form drag. This probably means that the velocity is increased 
in the wake for frequencies less than 1000 cycles per second. There seems 
little hope^ however^ of reducing flameholder drag by imposing high- 
frequency distiirbances of this kind. 

Turbulence intensity and spectra . - The turbulence component u’ 
in the direction of flow was measured with 0.0002- by 0.08-inch tungsten 
wires^ using the NACA constant-temperature anemometer amplifier. The 
wire was oriented parallel to the flameholder lips. No flame was present 
for these measurements. The results are presented in the following table: 


Nominal 
velocity, 
U, ft/sec 

Distance 

from 

nozzle 

lip, 

X, in. 

Distance from 
top wall, 
y> in. 

Vertical 

distance 

from 

station, 

in. 

Turbulence 
component, 
u*, ft/sec 

Scale of 
turbu- 
lence. 
Lx, in. 

50 

0 

1/2 ' 

Measured 

1/4 

0.306 


100 

0 

1/2 

J 

> from 
centerline 

1/4 

.420 


50 

6 

1/8^ 


1/4 

1.38 

0.2 

50 

6 

1/8 

Measured 

1/8 

1.68 


50 

6 

1/4 

from 

1/4 

.818 

.15 

50 

6 

1/4 

>flame- 

1/8 

.978 


50 

6 

1/2 

holder 

1/4 

1.021 

.15 

50 

100 

6 

6 

3/4 

l/8j 

lip 

1/4 

1/4 

3.35 

2.44 

.17 


The spectra measured at several stations were compared with spectra 
found representative for isotropic turbulence (refs. 35 and 36). A typ- 
ical plot is shown in figure 24. The scale is estimated as the one 

whose spectrum most nearly corresponds to the measured spectrum (refs. 

35 and 36). This scale^ where determined is given in the preceding 
table. 

The apparent transition from laminar to turbulent flow in the bound- 
ary layer on the flameholder surface occurred if the flameholder was 
about 1 inch or more from the test-section nozzle; the disturbance veloc- 
ity in the boundary was roughly that in the nozzle itself. With the 


^This is the scale that would be used for calculating a turbulent 
burning velocity from existing theory (ref. 8). Herein^ it is of interest 
because of the scale of disturbance the flame amplifies; ^ 1. 
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flameholder within 1 inch of the nozzle, eddies were found to shed at 
420 cycles per second when no flame was present. With the flameholder 
at 1 inch or more from the nozzle, transition occurred upstream of the 
flameholder lip and no clearly defined eddy-shedding frequency was found 
in the wake. 

With a flame anchored on the flameholder, no eddy shedding was ob- 
served without external excitation. The flame appeared turbulent for 
both flameholder positions greater or less than 1 inch from nozzle. It 
appears, therefore, that the transition to turbulence in the separated 
boundary affects the turbulent flame appearance in about the same manner 
as transition in the attached boundary layer does. The lean blowout 
fuel-air ratio was adversely affected when the transition was moved from 
the attached to the free boundary layer. 

In an effort to obtain a laminar flame zone, some exploratory work 
was done to move the transition point farther downstream of the point of 
separation from the flameholder. Small unheated cylindrical rods (as in 
ref. 3) placed about 1 inch from the nozzle entrance were used as flame- 
holders. It was possible to move the transition point about 1 inch down- 
stream of the flameholder (at 50 ft/sec). Lean blowout, however, was 
moved to a fuel-air ratio rich enough to cause excessive heating in the 
glass walls, which impaired the use of the shadowgraph system. In the 
transient period immediately after ignition while the walls were still 
cool, it was possible to observe a clear picture of the flame. In order 
to study the flame at length, this course was abandoned, and it was 
decided to accept the turbulent transition at or upstream of the lip of 
the flameholder. 


Calibration of Disturbance Amplitude Against Speaker Input 

Owing to the interaction between the standing wave tube supplying 
the sound to the flameholder and the standing waves in the test section, 
it was necessary to calibrate the sound input with the flameholder at 
the position and velocity to be used in the, burning tests. For most runs 
the flameholder was placed midway between acoustically open ends in the 
test section in order to suppress velocity disturbances at the flame- 
holder due to all even modes of vibration of the test section. 

The inputs were set in the following manner: The flameholder posi- 

tion and flow velocity to be used in a burning run were set. The hot 
wire was placed l/8 inch above the flameholder lip, l/8 inch upstream of 
flameholder lip so as to be centered over the port from which the sound 
issued. Resonant frequencies were excited in the flameholder support 
tube. The input to the speaker required to give the desired u* at the 
wire was recorded for each frequency. The wire was removed, fuel intro- 
duced and ignited, and the burning run was made. The distribution of 
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V* resulting from this sort of excitation was not explored. The change 
in the nature of interaction between test sections and the standing wave 
tube with and without combustion was not explored but was presumed to be 
small, since only about 30 percent of the combustible mixture was burned 
in the test section. 


St/Sl Determined by Photomultiplier Probe 

The author (ref. 37) and Kaskan (ref. 38) have used photomultiplier 
tubes to sense transient area changes in a laminar flame distorted by 
sound. Clark and Bittker (ref. 33) used photomultiplier tubes to study 
bunsen flames that were alternately laminar and turbulent by comparing 
the average intensity of light emitted from the flame. They found that 
turbulent and laminar flames having identical mixture ratios and feed 
rates produced identical quantities of light. Therefore, it appears 
that, if used with caution, the photomultiplier tube can become an 
average -heat -release -rate meter. 

In the following discussion, the local heat release rate is charac- 
terized by the ratio of turbulent to laminar flame speed Sij/Sj^. This 

ratio Sqi/Sl is defined as the heat released per unit time divided by 

the heat that a pair of laminar flames alined with the axis of the duct 
would release per unit length per unit time. The relation is illustrated 
by considering the way the quantities are measured. 



Region photo- 
multiplier scans 


Combustibles 


Sound 


Combustibles 
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In the test section shown in the preceding sketch, flame fronts are 
anchored on the top and bottom of the flameholder and appear to span 
from wall to wall. The photomultiplier scans a narrow region of the test 
section. The voltages read from the photomultiplier tube output are re- 
corded at the same fuel-air ratio for a turbulent flame and again for a 
laminar flame. The flame-front half angle is measured in each case and 
each voltage multiplied by the cosine of the half angle. The ratio of 
the resulting product for the turbulent flame to that for the laminar 
flame is defined as 8^1/8^. In the region of interest the cosine of the 

half angle was approximately 1 for the turbulent flames so that, in ef- 
fect, the unit length of the flame fronts becomes the unit length of the 
combustor . 

In obtaining the laminar flame, the same fuel and air rotameter 
settings were used as for the turbulent flame but a large fraction of the 
mixture was bled off upstream of the flameholder. Three banks of 200- 
mesh screen were installed just upstream of the flameholder and the re- 
sulting flame appeared laminar. Near the flameholder the flame did not 
span from wall to wall, so that the voltage recorded for this region was 
quite low. Far downstream from the flameholder the flame front became 
twisted so that it represented a flame depth greater than the wall-to- 
wall depth. In between these two regions was a region in which the flame 
appeared to touch the walls and have negligible twist. The voltage re- 
corded by the photomultiplier probe as it scanned this region was nearly 
constant. The extreme voltages recorded in this region envelop a band 
within which the true voltage corresponding to a pair of laminar flames 
was thought to lie. 

In order to pinpoint an exact value for the laminar voltage as well 
as to check the method (the zone of uncertainty is quite small), an aux- 
iliary determination of was made. Readings were made of 

photomultiplier-probe voltages at seven equally spaced stations along the 
combustor from flameholder to exhaust nozzle, and an average voltage was 
calculated. 8imultaneously, measurements of the momentum pressure drop 
across the flame zone were made. With the help of 8cur lock's results 
(ref. 1), the fraction of combustible burned F was determined from the 
measured momentum pressure drop. Then, from 

Length x 2 x 8m x depth 

F = i 

Uq X depth X height 


(where the length, depth, and height are the dimensions of the combustor) 
an average 81J was calculated. Using this 8ip, the measured average 

voltage, and the value for laminar flame speed 8^ taken from Dugger 
(ref. 39), voltage corresponding to a laminar flame was obtained from 
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Sl Srp 


where and Vip are the photomultiplier voltages for laminar and 

turbulent flames. This voltage lies within the zone of uncertainty found 
by direct viewing of the laminar flame. These results are reported more 
fully in the section Examination of the Method. 

It is- obvious that citing Clark and Bittker (ref. 34) does not imply 
that any photomultiplier calibration was obtained from their work. Ref- 
erence 34 established that (for a given combustor at a given composition 
and flow rate) a photomultiplier voltage is proportional to the heat re- 
leased within sight of the photomultiplier tube and is independent of 
the existence or lack of turbulence. The agreement of the two methods 
for determining seems to justify this use of the photomultiplier 

probe . 

It was hoped initially to follow the amplitude of a flame wrinkle 
by measuring the a-c component of the photomultiplier signal. The rea- 
soning was that the flame would be wrinkled sufficiently for its local 
instantaneous slope to serve as a detectible index, of flame-wrinkle am- 
plitude. This slope, of course, would show up as an increased voltage 
at the photomultiplier because the length of the flame cutting across 
the field of view would be increased. 

The data taken showed an amplitude more dependent on the distance 
from the duct exhaust than on the distance from the flameholder. It was 
discovered .that the sound waves introduced as the source of excitation 
at the flameholder reflected from the exhaust of the combustor. The a-c 
component read by the probe, therefore, was largely influenced by the* 
time -varying displacement of local flame zones due to plane waves of 
sound. Figure 25 shows a comparison of the filtered a-c signal for two 
frequencies as a function of distance from the flameholder. One set is 
obtained from the photomultiplier probe with burning, the other set from 
a hot-wire anemometer without burning. 

The use of the photomultiplier probe was restricted to surveying 
for the effect of dist\irbances on local S^/Sj^ values on a time-average 

basis only. 


RESULTS 

The analytical part of this paper suggests certain behavior that can 
be examined experimentally. Accordingly, tests are included that furnish 
the necessary data for a comparison with the results of the analysis. One 
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of the primary variables studied is the effect of the frequency of ex- 
citation on disturbance growth. The frequency range of interest happened 
to include some of the resonant frequencies of the combustor gas volume 
as well as the spontaneous eddy-shedding frequency of the flameholder in 
cold flow. Part of the tests studied the effect of combustor resonance 
or spontaneous eddy shedding on the results. 

Certain precautionary measures were taken in choosing geometry and 
operating conditions. To guard against burner resonance, the combustor 
was kept short so that only about 30 percent of the combustible was con- 
sumed. It was found that, when this value was exceeded, spontaneous 
burner resonance sometimes occurred. To guard against spontaneous eddy 
shedding, the conditions were kept within the Reynolds number reported 
in reference 1 to be free of spontaneous eddy shedding when a flame is 
anchored on the flameholder. Although these measures guarded against 
the occurrence of spontaneous oscillations due to the above two causes, 
they could not prevent the system from favoring certain frequencies with- 
in the range of interest. These influences, then, received special 
attention. 

Markstein has shown (ref. 14) that instabilities of a laminar flame 
front cause the front to assxmie a cellular structure under certain con- 
ditions. The occurrence of such cells clearly could confuse an inter- 
pretation of the results. Consequently, the data were taken solely at 
conditions where the cells do not spontaneously arise. 


Examination of the Method 

Runs were made at three fuel-air ratios at 50 feet per second and 
a range of disturbance frequencies. The flameholder position was 6 inches 
from the open end of the duct and 6 inches from the tunnel entrance. The 
0.306 -inch -wide gutter with antisymmetric disturbances was used. The 
change in impact pressure with and without combustion was measured by a 
single impact tube upstream of the flameholder. Since the difference is 
of the order of inches of water and the rate of mass flow is constant, it 
is assumed that the difference is entirely due to a change in static pres- 
sure. By use of Scur lock's results (ref. l) and the measured pressure 
drop, the fraction burned was calculated. Conciirrently, the light intens- 
ity at seven stations along the flame was measured and averaged. 

Figure 26 shows a comparison of the fraction burned with the average 
voltage (which is proportional to the average intensity) made by plotting 
both against the frequency of disturbance for the three fuel-air ratios 
examined. The following table demonstrates the amount of scatter in the 
ratio of fraction burned to mean voltage. 
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Fuel-air ratio 

Frequency, f, cps 

Average 

F/Vt 

0 

380 

470 

760 

1500 

F/Vt 

0.0444 

0.137 

0.131 

0.134 

0.137 

0.131 

0.134 

.0515 

.1017 

.1070 

.1065 

.1063 

.1035 

.1056 

.0580 

.0818 

.0817 

.0863 

.0858 

.0852 

.085 


From the relation 


p 2 Sqi X length x depth 
Uq X height x depth x 

(where lengthy height^ and depth are dimensions of the combustor) an 
average value of Sip/v^ is calculated. Using the average values in the 

preceding table for f/v^ at each fuel-air ratio^ the following values 
of Srjt/Vfji were obtained 


Fuel-air 

Srp/V qi. 

ratio 

(ft/ sec)/volt 

0.0444 

■ 1.115 

.0515 

.88 

.0580 

.709 


For these fuel-air ratios and the employed inlet temperature and pres- 
sure^ a laminar flame speed was obtained from Dugger (ref. 39). Then, 
from 


V 


^ Si^/Viji 


the voltage that would correspond to a laminar flame was calculated. 

In the following table are tabulated Dugger's flame speed values 
the value of calculated from the preceding relation, and the 

zone of measured directly by viewing a pair of laminar flames in 

the test section. 





NACA TN 3830 


55 


Fuel-air 

ratio 

Sl 

(ref. 39) 

Sl 

^ Sgi/Vip 

"by direct 
measurement 

0.0444 

0.74 

0. 663 

0.658 to 0.682 

.0515 

1.10 

1.25 

1.15 to 1.39 

.0580 

1.32 

1.862 



Sl 

The taken from values' of — that resulted from momentum- 

pressure -drop measurements falls within the limits of by direct 

measurement in the table and is used as the calibration of the photo- 
multiplier probe. Local values of Sqi/Sl reported subsequently are de- 
fined by the ratio of voltage read at that station to this Vl for the 
applicable fuel-air ratio. 


Comparison of Flame Speed with Flame Appearance 

Shadowgraphs were taken at the same conditions used in the preced- 
ing comparison between average voltage and Ap. Some of these shadow- 
graphs, along with the distribution of S^/S^ for each case, are pre- 
sented in figure 27 . 

In most of these photographs the gross disturbance caused by the 
input sound can be identified. If a line is drawn along the turbulent 
flame front faired through all but these gross -wrinkles, its length 
divided by the length of a line faired along the undisturbed surface 
should be equal to 

Area of flame disturbed at frequency f 
Area of undisturbed flame surface 

A similar ratio is available from the average voltages obtained by the 
photomultiplier survey in figure 26 defined by 




(Vt) disturbed 
(Vqi) undisturbed 


The following table shows a comparison of the (Ajj/Aq)^ obtained by 

the former method and (Ajj/Ay)p obtained from the photomultiplier 
readings . 
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Frequency, 
X, cps 

Fuel-air ratio 

0.0444 

0.0515 

\%/P 

\%/A 

\%/P 

( — 
\%/A 

0 

1 

1 

1 

1 

380 

.18 

1.15 

1.-09 

1.11 

570 

1.304 

1.185 

1.19 

1.12 

770 

1.23 

1.23 

1.10 

1.18 

1500 

1.10 

1.085 

1.09 

1.07 


In general the correlation is good. It may he significant that at 
570 cps the (Ap/Au)p is about 10 percent greater than (Ap/Ay)^. At this 

frequency^ it is relatively easy to measure the area; that is, experi- 
mental scatter should not be as great as 10 percent. The evidence sug- 
gests the generation of additional turbulence in the flame, over and 
above the growth of the low-frequency-imposed disturbance. It is regret- 
table that a laminar stabilized flame could not be obtained in this equip- 
ment, because this frequency could conceivably cause a transition point 
from laminar to turbulent flow. 

For each plot (fig. 27(a) to (e)) u* is a constant, so that u^ /^L 
decreases with increasing fuel-air ratio f/a. decreases, for the 

most .part, with increasing fuel-air ratio and, therefore, with decreasing 
u*/Sl. 


position, the ratio (Ap/Au)p defined by 

descends with ascending fuel-air ratio. This is in 

keeping with the interpretation in the STABILITY ANALYSIS. The following 
table gives the values of Ap/Ay at 570 cps for the three fuel-air 
ratios. 


At a given x 
( Srp/Sp^ ) ^ P s turb ed 
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in. 


Fuel -air ratio 



0.0444 

0.0515 

0.0580 


^D 

"^70 

^D 

^^^70 


^^70 



Sm 

■*•0 

% 

Sfp 

■‘■0 


Sm 

-^0 

1 

1.522 

2.025 

1.129 

1.220 

1.076 

1.055 

2 

1.183 

1.437 

1.205 

1.107 

1.122 

.962 

3 

1.166 

1.376 

1.175 

1.08 

1.119 

.937 

4 

1.177 

1.358 

1.151 

1.058 

1.147 

.913 

5 

1.197 

1.352 

1.129 

1.038 

1.116 

.873 


Distribution of Local S^/S^ Values with Respect 
to Various Frequencies 

Figure 28(a) shows a plot of Sqi/Sl against distance from the 
flameholder. Figure 28(b) is a cross plot of Sip/S^ against frequency 

at several stations. Great care was taken in keeping fuel-air ratio, 
velocity, and plate voltage to the photocell constant throughout the run. 
The 0.306-inch flameholder, modified to prevent leakage about the edges, 
was used. The disturbance velocity measured as described earlier gave a 
u* =2.16 feet per second for all frequencies. 

In order to compare these data with theory, the following computa- 
tions were made: 

(1) From the preceding calibration the fraction burned is calculated 
as a function of combustor length for the upper envelope of the Sip/S-^ 

curves for the disturbed flames. 

(2) An idealized '^ 2/^1 against length was obtained from this 
fraction burned using Scurlock's results (ref. l). 

(3) Increments of length over which a given frequency has a maximum 
St/Sl were read from figure 28(a) . 

(4) These increments were plotted against ^ 2/^1 compared with 

the theoretical plot taken from figure 11. Recall that 
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Results of steps (l) and (2) are shown in figure 29 , and steps (3) and 
(4) in figure 30 . 

Similar data were taken at the lower amplitude, u‘ = 1.44 feet per 
second, and gave similaj: results. 


Disturhance-Displacement Growth 

In this section the actual disturbance -displacement growth is com- 
pared with that of qualitative theory (fig. lO) . 

In an effort to produce low-disturbance amplitudes, the glass walls 
of the test section were reduced to 6 inches in length and the flame- 
holder moved to 5/8 inch downstream of the nozzle entrance. At this 
position the speaker inputs were calibrated With the hot wire as previ- 
ously discussed, and shadowgraphs were taken for several low amplitudes. 
As mentioned under "Turbulence intensity and spectra", the flameholder 
boundary layer underwent transition in the separated layer. The photo- 
graphs of figure 31(a) show that, at these low intensities, the apparent 
width of the flame was not appreciably altered. The mean flame width 
measured from a direct photograph was used to determine h 2 /h-j^. 

The displacement h was measured as the distance from surrounding 
crests to a given trough in the flame surface. Measurements were made 
on both sides of the flame for each photograph. No attempt was made to 
separate the component of h because of superposition of a random ve- 
locity disturbance upon the imposed disturbance. The data for two am- 
plitudes and two frequencies obtained in this manner are shown in figures 
31(b), (c), and (d) . The dashed curve is the qualitative curve from 
figure 10 . 

The knee should occur between (h^/h-) )^„y and (h2/h]_)2j.i^j^2a,l 
discussed in the section Relation of Theoretical Results to Measurable 
Quantities. The vertical dashed lines in figure 31 are drawn through 
these two values. It can be seen that, and with the general shape of 
the qualitative dashed curve known, if a line were faired through the 
data points, it would exhibit a knee in the region bounded by ■ (hp/h-| ) ma.y 

and (h2/h]_) critical* ^ comparison of figures 31(c) and (d) shows that 
the knee moves to increased values of h 2 /h 3 ^ as frequency is reduced. 
Both the value of h2/h]L where the knee occurs and the variation of 
this value with frequency agree with the results of the analysis. 


Disturbance Phase Velocities 


From the photographs of figure 31(a) , values of Cr/Up were deter- 
mined, again using the mean flame width to furnish h 2 /h]_. The quantity 
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Ui/Uq then determined from Tsien's results (eq. (3)). Wavelengths 

were measured from trough to trough, and the mean position of a pair of 
troughs located the applicable value of hg/h^^. The results for two fre- 
quencies are shown in figure 32 > 

The dashed lines are the values of Cp/Ui obtained from figure 5(b) , 
which is based on the results and interpretation of the stability anal- 
ysis. The scatter of the experimental data is quite bad, but the magni- 
tude and trends of the data remote from the flameholder are quite good. 
Values of less than unity are obtained for values of h2/h2_<1.5, 

which implies flow reversal within this region. 


Effect of Test-Section Depth 

A single run was made to determine the effect of the test-section 
depth on the results. Figure 33 shows the data obtained with the 1-inch- 
deep combustor at a velocity of 50 feet per second and a fuel -air ratio 
of 0.0493. The turbulence at the center of the duct was half that of the 
l/2-inch duct. Combustor depth seems to have some influence on the tur- 
bulent flame speed, but not on the effect of imposed disturbances. 

The undisturbed flame speed in the 1-inch duct arrived at almost 
twice laminar flame speed in a longer time than with the l/2-inch duct. 
The effect of imposed disturbance on amplification as a function of fre- 
quency remains the same, however. This can be seen from a comparison of 
figure 33 with figure 28(a) . The range over which a given frequency 
dominates, taken from this run, is also plotted in figure 30. 


Effect of Amplitude 


The effect of amplitude on single frequencies is shown in figure 
34 , where flame speeds for two levels of excitation (580 and 960 cps) 
are compared. For the lower frequency, the increase A(S^ /Stp ) varies 


directly as u* at the two conditions. For the higher frequency, the 
curves cross. Thus, for lower excitation the effect is felt farther 
along the flame than for higher excitation. 


This behavior is not data scatter and was checked repeatedly. The 
same sort of thing occurs at 1500 cps, although not so markedly. The 
reason for this anomalous behavior is not immediately apparent. 


Effect of Flameholder Size 


Figure 35 compares and flame appearance for two flameholder 

sizes at a velocity of 50 feet per second and a fuel-air ratio of 0.044. 



60 


NACA TW 3830 


Figure 35(a) shows data for a l/2-inch flameholder and a variety of ex- 
citation frequencies. The 580-cps excitation is again the dominant one, 
but not so markedly as before. Figure 35(b) shows the set of correspond- 
ing shadowgraphs. Figures 35(c), (d), and (e) compare values of Srp/S^ 

for 0.500-inch and 0.306-inch flameholders and give a shadowgraph picture 
for each case. The effect of flameholder size can be summarized as 
follows : 

(1) In each case, flame speed is noticeably decreased with the 
larger flameholder. 

(2) The impact of the 580-cps frequency is not nearly so great in 
the l/2-inch flameholder as in the 0.306-inch flameholder. 

The decrease in flame speed with increase in flameholder size might 
be explained by the theory that (l) the wake region is longer and (2) 
the growth of the higher frequencies is therefore somewhat inhibited by 
this long, relatively stable wake. 


Effect of Velocity 

Momentum-pressure-drop and flame-speed readings were taken over a 
range of frequencies, using a 0.306-inch flameholder, a velocity of 100 
feet per second, and fuel-air ratio of 0.0493. 

The results are shown in figure 36(a). The metering systems drifted 
slightly at this setting, making it difficult to hold fuel flow steady 
and causing some scatter in the calibration Sl. 

The most sensitive frequency at 100 feet per second was 1700 cps, 
which is about three times the frequency observed at 50 feet per second. 
Clearly the effect is not allied with a spontaneous eddy-shedding fre- 
quency. The fraction burned is not so great as at 50 feet per second. 
Also, the flame speed increased sharply up to the end of the duct. These 
two factors would tend to emphasize the higher frequencies within the 
limited duct length, and would seem to be directly related to the insta- 
bility of the flame zone. The single data point at 1700 cps is included 
in figure 30 . 

This run was directed toward finding a frequency other than 580 for 
maximum effect, mainly to prevent duct resonance. A different frequency 
for maximum effect was obtained. Shadowgraph photographs of these flames 
are shown in figure 36(b) . Again, the disturbance is antisymmetric. 

The possible effect of duct resonance on these results must be con- 
sidered. The modes of resonance likely to occur in a combustor are the 



NACA TN 3830 


61 


"organ-pipe" or longitudinal modes, the transverse modes, and a single 
mode coupling the combustion chamber with the volume upstream of the 
nozzle. In the configuration studied, the longitudinal modes have fre- 
quencies, roughly, of 250 or 500 cps plus multiples of 500. The trans- 
verse modes have a minimum frequency of 3000 cps, so they are not con- 
sidered here. The last type of oscillation gives a frequency of the 
order of 45 to 90 cps and is below the range of interest. 

A longitudinal mode could affect the results in two distinct ways: 
(l) If the flameholder were placed in a standing wave system at a pres- 
sure node, the time-varying velocity past the flameholder would cause 
large initial disturbances that would be symmetric with respect to the 
duct axis. (2) If the flameholder were placed at a pressure loop, there 
would be a coupling between the combustor and flameholder resonance that 
would be initially compensated for by the input calibration already dis- 
cussed. If the nature of duct resonance shifts the resonant frequencies 
so that this compensation is no longer adequate, the input u' is no 
longer constant for all the frequencies. 

Because the photographs of figure 36(b) show the disturbances to be 
antisymmetric, the first manner of coupling is rejected. To explore the 
second manner of coupling, the li5)ut u* for the frequencies on either 
side of the most sensitive ones (i.e., 1500 to 1700 cps; see fig. 36fal ) 
were increased by a factor of 3. In no case was it possible to increase 
the fraction burned as much as the 1500- to 1700-cps disturbances did. 

It is concluded that the results shown in figure 36(a) are due to the 
sensitivity of the flame zone to the imposed disturbance and are not due 
to duct resonance. 


Symmetric Disturbances 

A symmetric disturbance was used for tests with a 0.315-inch flame- 
holder, velocity of 50 feet per second, and fuel-air ratio of 0.044. No 
separate calibration of input sound was attempted, the Inputs of the 
antisymmetric-disturbance tests being used. 

The results are shown in figure 37(a) and corresponding photographs 
in figure 37(b) . Comparisons with flame speeds obtained with antisymmet- 
ric disturbances are given. There was a leak of approximately l/l6 inch 
past one edge, which provided an additional flame surface and increased 
the apparent mean flame speed near the flameholder but did not disturb 
the results at some distance from the flameholder. In general, symmetric 
and antisymmetric disturbances gave roughly comparable shifts in flame 
speed near the flameholder. At a distance from the flameholder, anti- 
symmetric disturbances increased the flame speed more than did symmetric 
disturbances. These results are in accordance with figure 10. 
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DISCUSSION 


This section is devoted to a discussion of the questions raised hy 
the analytical and experimental phases of the investigation. 

1. In vhat respects do the experimental results and the interpreta- 
tions'^ the stability analysis agree ? - The areas of agreement are sum- 
marized in the following table.: 


Interpretation of stability 
analysis 

Investigation of sensitivity of 
a flame in a duct to transverse 
velocity disturbance 

1. The flow is stable to symmetric 
disturbances, unstable to antisym- 
metric disturbances. A qualifi- 
cation is made near a flameholder 
where the two types of disturb- 
ances are equally unstable. 

1, The flame seems about equally 
unstable to symmetric and anti- 
symmetric disturbances near 
flameholder, is more unstable 
to antisymmetric disturbances 
far from flameholder. 

2. In moving along the flame at 
phase velocity, a disturbance of 
given frequency grows until the 
flame width becomes about 0.4 of 
the wavelength, then an apparent 
decay of the disturbance begins. 

2. A disturbance of given fre- 
quency first grows, then de- 
cays. The critical flame 
width where decay starts is 
slightly smaller than that 
predicted. 

3. There is a most -amplified fre- 
quency having a wave number 
2ctQh-L = 6 -*■ 8. 

3. There is a most -amplified fre- 
quency which, for the lengths 
Investigated, corresponds to 
2oQh^ “ 12. 

4. The phase velocity ranges from 
a velocity equal to the cold-gas 
velocity at the flameholder up to 
1.5 times the cold-gas velocity 
far from flameholder. 

4. The phase velocity behaves 
much as predicted except near 
the flameholder, where it is 
lower than the cold-gas 
velocity. 

5. The amplification depends on 
v’/^Tq ^ given dh2/dx. 

5. The amplification depends on 
v'/^Tq ^ given dhg/dx. 


The details of the system that would define more exactly the short- 
comings of the stability analysis as well as the adequacies await the 
development of a transient-velocity measuring technique that can be used 
in a flame zone. In at least one instance there was evidence that a form 
of transition took place (see Comparison of Flame’ Speed with Flame Ap- 
pearance). The stability analysis cannot indicate when this takes place. 
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2. Would the same distribution of displacement and Sip/Sj^ with dis- 

tance~from the flameholder obtain if no growth occurred in the flame 
zone ? - This seems unlikely. For example, if the input disturbance dis- 
tribution simulates a source at the corner of the flameholder, only near 
the flameholder would there be a marked transverse velocity component. 

The highest frequencies would show the greatest increases in flame speed. 

3. Do eddies shed ? - They certainly seem to. In figure 38 photo- 
graphs of a lean flame, both the shadowgraphs and the direct photographs 
taken through the stroboscope, are compared with a smoke-trace photograph 
of eddies shedding from a rod at a Reynolds number of 73 (ref. 32). The 
appearance of flame clearly indicates that eddy shedding is present. 

As the fuel-air ratio is increased, the range over which these eddies 
appear near the flameholder narrows. Figure 39 shows a plot of frequen- 
cies at which eddies appeared against fuel flow from lean to rich blow- 
out. These data were obtained by exciting the flame at frequencies at 
which the constant -speed stroboscope would permit observation, that is, 
approximately every 100 cps. 

Eddies are characterized by lobes of flame of about equal size on 
either side of the flameholder. These are contrasted to mere sinuosity 
of the flame. At 700 cps the lobes appeared symmetric, indicating duct 
resonance at this frequency. At 800 cps the lobes were antisymmetric 
and were observed throughout the fuel-air ratio range. At the lean fuel- 
air ratios, most of the disturbances produced eddy shedding. 

4. Does the turbulent -flame propagation velocity depend upon up- 
streair~turbulence? - It should. Bolz arid Burlage (ref. 11) and Mickelsen 
and Ernstein (ref. 12) have shown turbulent flame speeds to be in fair 
agreement with Sciir lock's theory (ref. 8), omitting flame -generated tur- 
bulence effects. These data were obtained in turbulent "soap-bubble" 
flames. Thus, in addition to the effect of the instabilities mentioned 
here, the incident turbulence should have an effect. 

The most important point, however, is that the instability merely 
acts as an amplifier for a relatively narrow range of frequencies, the 
output within limits of the linearized theory being proportional to the 
input. Part of the input is, of course, the turbulent motion. The other, 
and in some cases the most important part, is the role played by the 
sound in the duct which, upon hitting the flameholder, produces a trans- 
verse velocity that may be subsequently amplified. In that this sound is 
related to the incident tirrbulence, the turbulence has an effect. 
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Some data on an effect of incident turbulence were obtained while 
exploring the properties of the system. Figure 40 shows a plot of Sij/S-^ 

against distance from the flameholder x for two flameholder positions 
and two fuel -air ratios. The table in the section "Turbulence intensity 
and spectra" shows that the turbulence intensity increases from inlet 
nozzle to exhaust and from the wall to the interior of the test section. 
Taking as representative the value of u' measured 3/l6 inch above the 
flameholder lip, the values of the turbulence parameter u'/s^ vary 

over a factor of 2 for the two flameholder positions. It is apparent 
that, everywhere except near the flameholder, there is little effect of 
the initial u’/Sj^ values measured. The data eire insufficient for pin- 
pointing the cause of this lack of effect. Two surmises appear warranted: 

(1) The flame was noisier at the 9-inch position than at the 6-inch 
position because of flame-duct interaction, u* was measured in the 
absence of a flame. 

(2) The turbulence found upstream of the flameholder was due to 
transition in the boundary layer. The u* in the boundary layer was 
growing for both positions measured. It is thought that the growth of 
this u* is accelerated when the boundary layer separates from the 
flameholder and quite rapidly achieves a value that would be relatively 
insensitive to u’ prior to separation. 

Additional data on the effect of upstream turbulence were obtained 
by putting three banks of 200 -mesh screen just upstream of the flame- 
holder. Shadowgraphs of the flame in this case with and without a 560 
cps disturbance are shown in figure 41 . Also plotted in this figure is 
S^/Sj^ measured with and without* screens against distance x. In the 

absence of imposed disturbances, Srp/Sj^ is lower with the screens than 

without. Although the conditions were set for 50 feet per second, leaks 
were found between the mixing section and the flameholder for this run 
so that the velocity was not known. The reduction in Sgi/S-^ is thought 

primarily due to the reduction in u‘ by the presence of the screens. 

It is of interest to compare the flame-speed data obtained in this 
investigation with data taken in a similar apparatus. Figure 42 shows a 
comparison of Sij/Sj^ obtained by photomultiplier tube with ob- 

tained by Scurlock (ref. 8) by measuring a mean flame width. Because 
was taken as a property of the mixture, the values of Siji/Sj^ near 

the flameholder are low. The values of S^i/S^ obtained with the rela- 
tively high turbulence appear higher than those obtained by Scurlock. 

The possible displacement of the absolute magnitudes of due to 

the different methods used would not warrant interpreting these data to 
give a quantitative dependence on turbulence intensity. 
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5. How does Sip/S-^ obtained from area measurements compare with 
Srp / taken from photomultiplier data? - The shadowgraphs indicate con- 

siderable wrinkling of the flame surface. In addition, it appears that 
density gradients exist within the outer envelope of flame surface, im- 
plying heat release there. 

Since no single continuous flame surface is apparent, measuring 
flame area from these pictures seems impossible. For the sake of com- 
parison, however, some relative surface extension should be of interest 
to those who would formulate a theory of turbulent flame propagation. 
Such measurements were made in the following manner:. The negatives of 
a group of photographs were used to project an image of the flame on a 

screen. The image of a 6-inch test section measured 2^ feet in length. 

A map measure was used to find the length of the outer envelope of the 
group of flame images, ignoring the obvious density gradients within the 
outer envelope. Lengths were measured on the top and bottom surface of 
each flame, averaged, an d div ided by projected length. The resulting 
quotient is defined as (S^/Sj^)^. 

In the following table, the ratios so obtained are compared with 
the average Sij/S^ obtained from a photomultiplier survey (S^/Sj^)p. 


Fuel-air 

ratio 

Uq. 

ft/ sec 

Flame - 
holder 
size, 
in. 

Photo- 
graph 
from 
fig. - 

Excitation 
frequency, 
f, cps 

( 

( 

(St/Sl)p 

\ L/A 

\^L/P 


(St/Sl)a 

0.0444 

50 

0.306 

28(a) 

0 

1.436 

2 

1.392 

.0444 

50 

.306 

28(e) 

560 

1.77 

2.54 

1.430 

.0515 

50 

.306 

28(a) 

0 

1.40 

1.82 

1.30 

.0444 

50 

.5 

28(a) 

0 

1.36 

1.735 

1.28 

.0493 

100 

.306 

37(b) 

0 

1.65 

2.460 

1.490 


In general, the trends indicated by the (S^7^)p followed by the 

(S^/Sj^)j^. In view of the obvious internal structure and three-dimensional 

nature of the wrinkles, no argument for or against the wrinkled flame as 
an adequate turbulent flame model can be Introduced. 

6. How can flame speeds be predicted on the basis of these results ? - 
The most promising method of predicting flame speeds would be based on an 
incident and amplified spectra, as recommended by Markstein (ref. 14). 

For example: If the incident v' spectrum corresponds to an Lx of 

0.15 inch and each frequency is amplified as the theory states, the 
spectra at ^ 2 !^! - 0.25, for 1/K]_ « 20, would correspond to 
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the ones shown in figure 43. Between the two spectra would be a con- 
tinuously shifting spectrum. What is needed, then, is a model of flame 
propagation that would allow for a continuously shifting spectrum of 
turbulence. This model should have the advantage of allowing for the 
appearance of discrete frequencies that dominate the balance of the 
spectrum. 


An interesting point comes up in regard to spectra. It is seen ex- 
perimentally as well as theoretically that very little effect is derived 


from exciting the flame at a frequency less than 


2 _U_ 

3 2h^ 


or greater than 


4u/2h3^ (see figs. 26, 36, and 7) . It would appear, therefore, that for 
the present conditions (i.e., velocity of 50 ft/sec, 2-in. duct) the 
energy in the spectrum of turbulence below 200 cps will not increase 
flame spreading. 


A plot of the percentage of energy in the turbulence at all frequen- 
cies below f against fL^/U from Dryden (ref. 40 ) is shown in figure 

44 . For an assumed scale Lx = frequency below the frequency 


2 U represents about 35 percent of the total turbulence energy. This 

3 2h^ 


amount of energy is contained in disturbances whose frequency is too low 
to be amplified to any usable extent by the flame zone. Similarly from 

4U 

figure 44, the energy in the spectrum below the frequency f = 


^or 


U 


2h^ 

« 0.66) is 85 percent of the total turbulence energy. Thus, 




about 15 percent of the turbulence energy is at frequencies too high to 
be appreciably amplified by the flame zone. 


The scale of isotropic turbulence that would put the largest 

1 


2 U 4U 

centage of energy in the frequency range ^ < f < can be 


from the extremum of the quantity (ref. 38) 


per- 

found 


tan“^ 


^ 4U Lx 1 


- tan“^ 






or 
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In this case the amount of energy within the frequency range is 55 per- 
cent of the total as compared to 50 percent for Lj^ = ^ 2h-j^; 40 percent 
for ^ 2h]^; and 45 percent for ^ ^ 2h-jL* 

7. Can the results concerning frequency be applied to multiple 
flameholders ? - If n flameholders are arranged so that the flow field 
can be divided into n identical segments, Scurlock's results hold for 
each segment. The results in the theoretical sections of this report 
for the stability of the flow field apply to each segment and rely on 
the condition that v’ vanishes at the boundary of each segment. Al- 
though one of several modes of instability may satisfy this condition, 
it is obvious that the foregoing results do not give any information 
about the other modes. 

To illustrate the occurrence of multiple modes, consider a profile 
examined by Lin (ref. 24, pp. 219-220) where 

U = A + B sin y 

Lin shows that the wave numbers for neutral stability of the different 
modes available are dependent upon the number of inflection points in 
the flow and can be given by 



The choice of Z = Xy yields the following profile which is similar 
to the one studied earlier in this report for a single flameholder 



For this case neutral stability occurs at aXy = aZ = jt^/3, 0. 
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If two flameholders are present, the profile (following sketch) may 
be approximated by setting I = 2Xy, 



For this case the wave numbers for neutral stability of the foiir modes 
obtained are 

^ I pr 'JT r. 

aXy = a 2 = « 0 

Two new modes appear in addition to the two for the single profile. 
Similarly^ for three flameholders^ 



that is^ there appear four roots in addition to the two for I = ly. 

In short, the single mode investigated in this report is applicable 
to the case of multiple flameholders, but only as one of several modes 
that combine to give the entire instability. 

8. Can the excitation at the flameholder permit increased combustion 
rates without destroying the burner walls ? - Yes. In these tests the 
sound level at the exhaust plane was only 93 decibels with the flame ex- 
cited. When the flame is excited by the destructive duct resonance, the 
sound level is of the order of 200 decibels (ref. 19). 

In order to apply this type of excitation to an engine, the control 
exerted herein (a speaker system) would not be needed. Instead, the dis- 
turbance can be introduced by suitably placed whistles, mounted on the 
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flameholder and driven by the flowing stream. The whistle ideally should 
have a number of strong overtones and a fundamental equal to the most 
sensitive frequency UQ/Sh. The overtones would serve to increase the 

heat -release rate near the flameholder^ while the fundamental would in- 
crease the heat-release rate far from the flameholder. 

9. What would be the sacrifice in stability and pressure drop for 
such an excitation ? - Both these penalties require further investigation^ 
but from work that has been done the price does not seem high. For ex- 
ample: When the flameholder was placed S/S inch from the nozzle entrance 

and a 6-inch-long test section was installed^ lean blow-out at 50 feet 
per second occurred at a fuel-air ratio of 0.042. At a fuel-air ratio of 
0.0453^ the disturbance amplitudes required to cause blowout were found 
for several frequencies. The results are shown in figure 45 . Minimum 
disturbance velocity was required at 780 cps. At this condition the 
flameholder without flame had an eddy-shedding frequency of 420 cps. The 
Strouhal number is thus 0.206^ agreeing well with the results of Roshko 
(ref. 36). The cold-duct resonance fundamental was 1000 cps. The rea- 
son for this result is not known. The main point here^ however^ is to 
show the effect of disturbances to be used for promoting flame propaga- 
tion upon flameholder stability. It was found that^ at a fuel-air ratio 
of 0.047^ it was impossible to induce blowout with speaker inputs as much 
as 10 times those used in the earlier tests. It seems safe to say that 
the imposed disturbances will narrow the stability range slightly but not 
prohibitively and the amount of narrowing will depend upon the frequen- 
cies imposed. 

An interesting observation was made in the study of lean blowout. 

As the lean blowout fuel-air ratio was approached^ the flame emitted a 
low-frequency noise at 87 cps. The average heat release rate dropped 
abruptly. A survey taken with the photomultiplier probe showed the flame 
to be intermittently extinguished for all positions greater than 7/S inch 
downstream of the flameholder. Typical wave forms observed at several 
stations downstream of the flameholder are shown in the following table. 


Distance downstream 

Wave form of photo- 

of flameholder^ in. 

multiplier signal 

7 

8 


2 

4 

v/\/v. 

vAV 

\JvAJ 
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Random photographs taken at this condition show intervals during which 
no flame was present in part of the combustor. These photographs appear 
in figure 46 . 

It is curious that the flame surface facing upstream at the break 
in the flame is strong and well-defined both in the photomultiplier and 
the shadowgraph observations. The flame surface facing downstream is 
conversely ill-defined in both methods of observation. Near the point 
of rupture the photomultiplier signal is approximately symmetric. No 
explanation for this is apparent. The most significant feature of these 
data is that this form of resonance accompanied by intermittent blowoff 
downstream of the sheltered zone was the only form of excitation that 
reduced the heat-release rate in the combustor. On this basis it is 
tempting to generalize that all oscillations that do not cause flame 
parture must increase the heat release rate in the combustor. 

For the frequency of maximum amplification (560 cps)^ an increase 
of as much as 30 percent in heat release or about 0.2 percent of the 
kinetic energy of the flowing stream KEstream obtained in the 6- 

inch-long combustor for an energy flux 

^sound _ 6 

^stream 2500 

In this case the noise level at the exhaust was 93 decibels. 

Another way of increasing the heat release rate that is indirectly 
suggested by the present paper is to locally accelerate the flow or alter 
the profile to a more unstable one by use of a convergent -divergent com- 
bustion chamber. One way these can be combined with the introduction 
of a V* at the desired frequency is to construct a combustor with a 
wavy wall. Using this method of excitation would require a new experi- 
mental determination of the desired wavelength and amplitude of the 
wall’s waves ^ but would eliminate any effect on flame anchoring. 

10. Can anything be done about predicting the behavior of axially 
symmetric flames in cylindrical ducts ? - A little can be done. Firsts 
the steady-state problem, solved by Scurlock, of the cylindrical flame 
in a cylindrical duct yielded a profile as opposed to the _A_ 

profile for plane flow. If the arc part of the curved profile is assumed 
parabolic, the cylindrically symmetric disturbance equation can be solved 
by elementary functions for the case of symmetric disturbances. This, 
however, accomplishes nothing, because the profile is neutrally stable 
to symmetric disturbances. This is the same result obtained for the 
_A_ profile in the plane case. For the case .of antisymmetric disturb- 
ances, the problem is more complicated and no simple solution is apparent. 
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If the profile is affected hy a flameholder in such a way that it 
may more neaxly be approximated by a _/ — profile composed of three 
straight-line segments and two parabolic arcs, a range of instability 
for symmetric disturbances is obtained. This profile can be compared 
with the. corresponding profile in plane flow Such an analysis 

is given in the appendix. A few calculations are given to permit com- 
parison with the data of figure 15, for conditions that h]_ approaches 
infinity and p is constant throughout. The results are compared in 
the following table: 




1 b* 

2ab 

2 Id 

Plane case 

Cylindrical 

case 

0.1 

2.26 

2.08 

.2 

2.55 

2.42 

1.0 

2.55 

2.54 


In these cases the critical wave numbers agree quite closely. At a 
first guess it would seem that the results for the case of plane flow 
would be a fair approximation for the case of cylindrical flow. 

In order to obtain a rough experimental check on this supposition, 

a li-inch-I.D. glass tube and its accommodating nozzle were installed in 

place of the rectangular test section. The l/2- by 0.306-inch antisym- 
metric flameholder was retained. At 50 feet per second and a fuel-air 
ratio of 0.047, resonance-free operation was obtained at a combustor 

length of 4 :^ inches or less. With the flameholder set to give this 

length, the disturbance inputs were calibrated as before to give a dis- 
turbance velocity of 1.8 feet per second. Then the burning run was made 
with the settings indicated by calibration. The flame speeds were meas- 
ured by the photomultiplier probe plus a measurement of an average local 
flame diameter. The results are shown in the following table: 
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Frequency^ 
f, cps 

^t/^l 

1 In. from 
flameholder 

2 In. from 
flameholder 

4 In. from 
flameholder 

0 

2,27 

3.31 

3.58 

200 

— 

3.56 

3.77 

385 

2.59 

3.56 

3.95 

580 

2.68 

3.56 

3.95 

770 

2.82 

3.65 

4.25 

980 

2.89 

3.70 

3.80 

1170 

2.89 

3.70 

3.73 

1360 

2.96 

3.65 

3.73 

1740 

3.09 

3.65 

3.99 


The flame speed in this case is almost twice as high as in the plane 
turner. Part of this increase is undoubtedly due to the disturbances in- 
troduced by the corners of the flameholder. The fraction burned in this 
case is 20 percent at the position 4 inches from the flameholder^ corre- 
sponding to h^/h^ = 0 . 7 • Therefore^ the frequencies in 4 inches of 

this duct may be compared with those obtained in 6 inches or more of the 
rectangular duct. The values of giving the greatest increase in 

flame speed at the termini of the respective ducts are 


Frequency^ 
f, cps 

2hi, 

in. 

Uq. 

ft/ sec 

Geometry 

2(XQhj^ 

for max. 
®t/^L 

570 

2 

50 

Plane 

11.95 

770 

1.5 

50 

: 1 

Cylindrical 

12.10 


Of course^ these data were obtained in a single test^ performed to 
satisfy curiosity. The agreement does make it seem likely that a 
carry-over of about 1:1 would be feasible in going from the plane case 
to the cylindrically symmetric case. 


SUMMARY OF RESULTS 
THEORY 

The flow field of a flame in a duct is unstable to transverse dis- 
turbances having wavelengths 2.5 times the local flame width. 
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The amount a given frequency disturbance will grow varies directly 
as the frequency and inversely as the rate of flame spreading dh/dx 
and is independent of velocity and density ratio across the flame front. 


f = 


The frequency where maximum velocity-distiarbance growth obtains is 


2.8 ^0 
rt 2hi^ 


where 


Uo 


is approach velocity and 2h-|_ 


is duct width. 


This maximum growth occurs by the time the flame fills 64 percent 
of the duct and implies an amplification v’/vq = (1.108)^/^1 where 
= dh/dx. Frequencies higher or lower than this are amplified less. 
Higher frequencies achieve their terminal amplification at smaller val- 


ues of h 2 /h]^. For example, a frequency of 

tion of v'/vq = (l.07)^/^l, where K-^ = dhg/dx 
fills 28 percent of the duct. 


has an amplifica- 
by the time the flame 


By using a crude model of flame propagation, it is possible to pre- 
dict the qualitative effect that amplified disturbances of discrete 
frequencies will have upon turbulent flame speed. 


EXPERIMENT 

The effects of disturbances imposed at the flameholder are in sub- 
stantial agreement with those predicted by theory. In addition, the 
results show: 

1. Eddies shed from the flameholder with flame when suitably 
excited. 


2. The turbulent- to laminar-flame-speed ratio varies inversely as 
flameholder size, directly as ratio of flow velocity to flame speed or 
ratio of disturbance velocity to flame speed. 

3. Potentially, the most useful theory of turbulent flame propaga- 
tion would be based on an incident and locally shifting spectrimi of 
turbulence . 

4. Gains in local heat release rates can be obtained by using acous- 
tic disturbances at a flameholder without damage to walls or excessive 
energy expenditure. 


CONCLUDING REMARKS 


The form of instability that is discussed in this paper is of in- 
terest in formulating a theory of turbulent flame propagation along with 
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plane combustion wave instability, and the boundary layer instability, 
mentioned earlier. At present it is not known whether the latter two 
forms of instability can be encouraged to influence turbulent flame 
propagation more than they do spontaneously. The little evidence that 
is available makes it appear that they cannot. The type of instability 
that is discussed in the present paper is of additional significance for 
the reason that its effect on the ratio of turbulent to laminar flame 
speed can be increased in some important applications. Several methods 
of increasing this ratio are as follows: 

1. Increase the disturbance velocity at the frequency the flame zone 
can amplify. 

2. Increase the amplification rate by locally raising the mean ve- 
locity in the combustor. 

3. Increase the amplification rate by locally varying the velocity 
profile. 

A method of employing all three of these improving factors would 
be to build a combustor having a wavy wall. 

Although the second and third methods listed will clearly cause an 
increase in pressure drop, the first method substantially increased the 
heat release rate at a cost of a very small fraction of the kinetic 
energy in the approach stream flow and a slight reduction in the operat- 
ing range of the flameholder. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, July 23, 1956 
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AFPEMDIX - AXIALLY SYMMETRIC DISTURBANCES IK 
AXIALLY SYMMETRIC FLOW WITH DENSITY JUMPS 

This appendix examines the stability of the flow field arising from 
a flame in a duct for the case of axial symmetry. The results are \ised 
in the body of the text to compare wave numbers for neutral stability of 
comparable plane and axially symmetric flow fields. The symbols are 
defined where tised. 


DIOTURBANCE EQUATION 

Pillow (ref. 25) defines the quantity as 

= <p(r)e^°^(2-ct) 

The velocity in the z-direction is 

u = U + u’ 


In the r-direction. 


V = v' 


(Al) 


where 


. V’ 



u' = 1 Mdii 

r dr 


(A2) 


and the continuity equation is satisfied identically, 
the disturbance equation then becomes 


where 



a^) q> - <p 




r dr 


1 

r2 


Pillow shows that 
i = 0 (A3) 


(A4) 


and R is Reynolds number. 

If the study is restricted to cases where R -► *, and U" = — U' 

^ r 

(i.e., parabolic profiles or flat profiles), then equation (A3) becomes 

(M - a2)<p = 0 
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or 


with solutions 


d^q) 1. ^ 

^2 r dr 


1 ? rx 

-?r + = 0 


(p = aj_l2(or) + a^K-^Car) 


(A5) 


(A6) 


where and are modified Bessel functions of the first and second 

kind. Therefore, 


V-' 


Silf 

= 5^ = 


iacp(r)ei°^(2-ct) 


= _ioei°'(z-ct) + agK^Car)] 


(A7) 


and 


u' = 


1 

r or 


1 . 

or r 


= (pt + iipeio-Cz-ct) 


= ei^-Cz-ct) j^al^(ar) + p l3^(or)j + ag j^aK]^(ar) + i K2^(or)j|. 


(A8) 


where the prime denotes differentiation with respect to the argument. 


BOUNDARY CONDITIONS 


Consider the following profile and density distribution 
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Let solutions to equation (A5) be denoted by <P 2 , <P 3 for regions 1, 

2, and 3 with densities pg, and P 3 , respectively. 


As in the plane case, the requirements are 


and because of symmetry 


Similarly, 


Since, at r = hj^. 


•e 

H 

11 

0 

at 

hi 

<Pl = <p2 

at 


(pg = <P 3 

at 

to 

(P3 = 0 

at 

r = 


" w 

I ) 

\^^/2 \^/ 


at he 




at h-: 




(A9) 


(AlO) 


^ ^ 2^.1 ^ 
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ai = -ai 


Il(alii) 
Kx (oiii) 


(All) 


At r = h2, <Px = <P2 with equation (All) 


<Pl 


[l(ohi) ^ 

Cx'Ccch^^) 


= (p 2 = 1)2^12^(0112) + h 2 Kx(ali 2 ) 


(A12) 


At r = hg, cpg = (Pj or 


(P 2 = b2^l2^(ah3) + b2Kx(ah3) 

= (P3 = c^I^(oh3) + C2K2_(oh3) 


(A13) 


Since K 2 _(otr) -»■ « as or 0, l 2 ^(or) -*-0 as or -*• 0, and at r = 0, 

<P3 = 0, 


C2 = 0 


(AM) 


Since 


^ = "p vsr ^ ^ 


' , du\ 

^ dr/ 


the pressure conditions yield at h 2 (from eqs. (A7 ) . and (A8)) 


-p, (-OC 


+iaU^) (acp* + ^ = -pg [ (-iocwu^) (cxcp^ + j| §] 


(A15) 


and at 




(A16) 
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BOUNDARY VALUE PROBLEM 
Employing the following notation, 

(oh^) 

ai = IiCohg) - KiCohg) 

Pi r 1 ^1 (°^i ) r 

■"2 = - K]15H5T 

p2 = K^CoLg) 

P3 Il(<^3) 


P4 — K]_(<^3) 


3e - ^iC^^eUdr/h 




Bq = Ki(ah2)(^ 


h2 


P9 = ^(al-)(ah3) ^ 

PlO = — (aKj^)(oh3) + ^ K3^(oh3) 
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and 


Ti = l^iah^) = 33 

Tg = al* (0213) + ^ 1^(0113) 


setting = 0 (this merely alters the real part of c hy U^) 
the determinantal equation for a^, bg, and is 


L 


ai 

0 

-cog 

0 


-Pi 

H 

cPs + Pe 
(113-0)39-321 


-P2 

P4 

cP? + Pg 

(U3-c)P2o-3ii 


^1 

0 

-ri 

0 

-(U3-C )rg 


= 0 


(AI8) 


One immediate result of equation (A18) is to show that if dU/dr 

vanishes at either h„ or h_ in region 2, c is real for all a'. 

u o s 

This can be seen by setting 


Pe = P8 ^ (^h2 = ° 


and factoring c from the third row of equation (A18), or setting 

Pll - H2 = = 0 

and factoring 113=0 from the fourth row of equation (A18). Either 
operation leaves a determinantal equation in which c is wholly real. 


If 



(A19) 


and 
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§2^2 = ^8 




^2 \ 
^ " ^11 ' 


PZ 

— y 


then equation (A18) can he reduced to 


co^ - 

ai 


Pi + cPg 


cci2 - g20'i 


an 


Pg + cPy 


(U3-c)t 2+— g3n 


ri 


P3+(^3"^)P9 ■ 


(U3-c)r2+— gsn 


n 


P4+(U3-c)pio 


This is in the form 


me + n 
qc + s 


oc + p 
tc + u 


= 0 


or 


mte^ - oqc^ + c(mu + nt - os - qp) + nu - ps 
and the discriminant of c is 


= 0 


D = m^u^ + n^t^ + o^s^ + q2p2 + _2munt - 2osqp - 
2muos - 2muqp - 2ntos - 2ntqp + 4oqnu + 4mtps 


(A20) 


= 0 


(A21) 


(A22) 
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If D < 0, then c is complex and an unstable mode exists, 
symbols of equation (A22) are defined 


The 


m 


Og 

P5 - Pi 57 




n = Pig 2 

o = P? - P2 5:^ 

P = P 2§2 


'I = P 3 n ■ ^9 


P9 T9 

s = U3P9 ■ ^ S3P3 - U3 — P3 


t = P 4 ^ - 3l0 


P2 ^2 

u = 0383^0 " ^ S 3 P 4 “ U 3 — 34 




(A23) 
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Figure 1. - Geometry of flow field. U, flow velocity in x-directionj p, densi 
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(b) Density ratio 5. 

Figure 4. - Concluded. Stability map for flame in a duct. Contours of constant amplification 

2b 

rate oc^ are given by lines of constant aCj^ Lines of constant frequency show path 

given disturbance originating at flamebolder takes as it propagates along flame zone at phase 
velocity c^.. 
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Figure 5. - Phase-velocity parameters. 
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Figure 7 . - Gains for flame filling the duct 





Figure 7. - Concluded. Gains for flame filling the duct 
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Figure 9. - Disturbance- velocity distributions. 
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Figure 15. - Effect of velocity profile on air^lif ication-rate parameter ac^ and on 
wave numbers for neutral stability and maximum instability. 
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Figure 20. - Schematic drawing of shadowgraph apparatus. 
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Figure 21. - Time-average velocity profiles with cold flow. 
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Figure 22. - Velocity profiles In wake of 0.306-inch flameholder 
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Figure 24. - Turbulence u' spectra for flow velocity U of 50 feet per second. 
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measured l/8 inch above flameholder lip. 
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Figure 25.,- Comparison of hot-wire and photomultiplier 
measurements showing effect of plane waves near exhaust 
of duct. Flow velocity U, 50 feet per second. 


Fraction 'burned 




(c) Fuel-air ratio, 0.0444. 


Figure 26. - CoiEparison of fraction 
burned obtained by momentum-pressure- 
drop measurement with average voltage 
obtained by photomultiplier survey. 

Flow velocity U, 50 feet per second; 
length, 6 inches; disturbance velocity 
u’, 2.38 feet per second; antisymmetric 
disturbance . 


Average volts 


Turbulent- to laminar -flame -speed ratio, S^/S^ 
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Distance from flameholder, x, in. 

(a) No disturbance. 

Figure 27. - Comparison of S>j/Sl and flame appearance 
for three fuel-air ratios. Flow velocity at plane of 
flameholder U , 50 feet per second; 0.306-lnch flameholder. 
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(b) With antisymmetric dlsturbeuice at 380 cps. 

Figure 27. - Continued. Comparison of Sq^/Sl and flame 
appearance for three fuel-air ratios. Flow velocity at 
plane of flameholder Uq, 50 feet per second; 0.306 -inch 

flameholder. 
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(c) With antisymmetric disturbance at 570 cps. 

Figure 27. - Continued. Conqparison of S^/S^ and flame 
appearance for three fuel-air ratios. Flow velocity at 
plane of flameholder 50 feet per second; 0.306-inch 


f lameholder . 
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(d) With antisymmetric disturbance at 760 cps. 

Figure 27. - Continued. Con?)arison of S^/Sl and flame 
appearance for three fuel-air ratios. Flow velocity at 
plane of flameholder Uq, 50 feet per second; 0.306-inch 


flameholder. 
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(e) With antisymmetric disturbance at 1500 cps. 

Figure 27. - Concluded. Comparison of Sqi/Sl and flame 
appearance for three fuel-air ratios. Flow velocity at 
plane of flameholder Uq, 50 feet per second; 0.306- inch 


flameholder. 
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(a) Variation of S^i/S^ with distance from flameholder. 


Figure 28. - Effect of frequency on flame-speed distribu- 
tion. Flow velocity U, 50 feet per second^ fuel-air 
ratio^ 0.0444; disturbance velocity, 2.16 feet per second 
test section depth, l/2 inch; 0.306 -inch flameholder; 
antisymmetric disturbance. 
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(■b) Variation of with frequency 

for various distances from f lameholder • 

Figure 28. - Concluded. Effect of 

frequency on flame- speed distribution. 
Flow velocity U, 50 feet per second; 
fuel-air ratio, 0.0444; disturbance 
velocity, 2.16 feet per second; test 
section depth 1^ inches; 0.306-inch 

f lameholder; antisymmetric disturbance. 
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(b) lig/^l agQ’i^ist X. 

Figure 29 . - Distribution of fraction burned F and idealized 
hg/hi for data of figure 28(a) . 
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Figure 30. - Frequency giving maximum increase in flame speed as function of hg/^l 
compared with results of stability analysis. 
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(a) Flame shadowgraphs for disturbance velocity u* of 0.75 foot 
per second at various frequencies f. 

Figure 31. - Displacement measured on flame front for small- 
amplitude disturbances. Flow velocity U, 50 feet per second; 
fuel-air ratio, 0.0484; turbulent disturbance velocity u*, 
0.303 foot per second. 
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(b) Disturbance velocity u', (d) Disturbance velocity u’, 

0.6 foot per second; fre- 0.75 foot per second; fre- 
quency f, 725 cps. quency f, 545 cps. 


Figure 31. - Concluded. Displacement measured on flame front for small- amplitude 
disturbances. Flow velocity U, 50 feet per second; fuel-air ratio 0.0484; 
turbulent disturbance velocity u*, 0.303 foot per second. 
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Figure 33. - Fleune-speed distribution in 1-inch-deep test 
section. Flow velocity at plane of flameholder Uq, 50 
feet per second; fuel-air ratio, 0.0493; l/4-inch flame- 
holder . 
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Figure 34. - Effect of aniplitude on flame- speed distribu- 
tion. Flow velocity U, 50 feet per second; fuel-air 
ratio, 0.0444. 
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(a) 0.500-Inch flameholder with antisymmetric disturbances. Disturbance 
velocities: 2.22 feet per second) Uturbulent^ 

second; iifpequency^ P®^ second. 

Figure 35. - Comparison of S^/S^ and flame appearance for 0.500- and 
0.306-inch f lameholders . Flow velocity U, 50 feet per second; fuel- 
air ratio, 0.0444. 
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Distance from flameholder, x, in. 

(b) Both flameholders without disturbance. 

Figure 35. - Continued. Coii 5 )arison of Sip/Sj^ and flame 
appearance for 0.500- and 0.306-inch flameholders. Flow 
velocity U, 50 feet per second; fuel-air ratio, 0.0444. 
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(c) Both flameholders with disturbance at 570 cps. 

Figure 35. - Continued. Conqparison of Sip/S^ and flame 
appearance for 0.500- and 0.306-inch flameholders. Flow 
velocity U, 50 feet per second; fuel-air ratio, 0.0444. 
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Distance from flameholder, x, in. 


(d) Both flameholder 8 with disturbance at 
1500 cps. 

Figure 35. - Concluded. Comparison of S^/Sj^ and flame 
appearance for 0.500- and 0.306-inch f lameholders . Flow 
velocity U, 50 feet per second; fuel-air ratio, 0.0444. 
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Distance from flameholder, x, in. 

(a) S»n/&L distribution and comparison of fraction burned 
with average voltage. 

Figure 36. - Sij»/Sl distribution, shadowgraphs, and com- 
parison of fraction burned with average voltage. Flow 
velocity U, 100 feet per second; fuel-air ratio, 0.0493 
antisymmetric disturbance. 
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Distance from f lameholder, x, in. 


(a) Sj/S^ distribution for symmetric disturb- 
ances. Disturbance amplitude unknown; input 
same as for figure 28(a) . 

Figure 37. - Comparison of obtained 

with flamebolder for producing symmetric dis- 
turbances with S^/S^ obtained with flame - 

holder for producing antisymmetric disturb- 
ances. Flow velocity U, 50 feet per second; 
fuel-air ratio, 0.0444. 
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(b) Compaurison of symmetric and antisymmetric flameholders 
with and without disturbance. 

Figure 37. - Concluded. Conqparison of Srp/S^ obtained 

with flameholder for producing symmetric disturbances 
with obtained with flameholder for producing 

antisymmetric disturbances. Flow velocity U, 50 feet 
per second; fuel-air ratio, 0.0444. 
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Shadowgrapli of excited flame; frequency f, J60 cps 


C-41538 


Eday trail 


Figure 38. - Comparison of photographs of excited flame with 
eddy trail. 




Relative fuel flow 







0 123 456789 

Distance from flameholder, in. 


("b) Fuel-air ratio, 0.0515. 

Figure 40. - Comparison of relative flame speeds at two flameholder positions. 
Flow velocity U, 50 feet per second. 
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Figure 41. - Con 5 )arison of flame speeds with suid without l)anhs of 200- 
mesh screens upstream of flameholder. 
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Figure 42. - Comparison of S^/S^ taken in 1/2 -inch- and 1-inch-wide ducts with 
Scurlock*s data. 
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Figure 43. - Spectra of incident and flame an^jllfied turbulence. 
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Figure 44, - Turbulence energy distribution for exponential correlation coefficient. 
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Figure 46. - Spontaneous oscillations showing intermittent 
flame interruption. Frequency f , 85 cycles per second; 
fuel-air ratio, 0.0444; flow velocity U, 50 feet per second. 
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